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ABSTRACT 



Specific heat measurements on three samples of polyethylene, 
differing only in density, were made in the temperature range of 2.U 
to 30°K. A definite density dependence was noted for the specific 
heat in this temperature interval which allowed extrapolation of the 
data to completely crystalline and completely amorphous cases. At 
the lowest temperatures, the amorphous results were observed to dis- 
play an "excess" heat capacity which could be accounted for by the 
occurance of a single delta function peak in the low frequency part 
of the vibrational spectrum. This excess did not appear in the com* 
pletely crystalline extrapolated data, and the specific heat was found 
to be proportional to the cube of the temperature up to 9°K. 

An attempt was made to compare the results with two theoretical 
models (Tarasov and Stockmayer-Hecht) with only fair agreement in one 
case, and none in the other. The agreement with previous experimental 
results is excellent if the density dependence and the excess heat 



capacity are considered. 
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1. INTRODUCTION 



General Remarks 

The objective of the investigation which is the subject of this 
thesis was to determine the specific heat of polyethylene by direct 
calorimetric methods in the temperature range between 2.4 and 30°K. 
This study was prompted by some earlier results of low precision which 
indicated a rather sizeable density dependence of the specific heat in 
the low temperature range [34] . This observation was in contradiction 
to statements in the literature which denied such a dependence [52j . 

In addition, theoretical models had predicted an unusual behavior of 
the specific heat of fibrous solids, such as polyethylene, at low 
temperatures which had not been verified by previous measurements 
[19, 46, 47, 52 ] . It was hoped that these effects, if they were to 
occur, would be found in the range covered by the present experiment 
thus affording a comparison of the experimental results with the pre- 
dictions. Further, it was desired to investigate the density depend- 
ence over a larger temperature range, and with more precision than 
that in the work noted above. 

The unusual heat capacity behavior which had been predicted to 
take place is a consequence of the polyethylene molecules being long 
chains, and the forces along the chains being touch larger than those 
between chains. This one dimensional characteristic gives rise to 
theoretical descriptions for the specific heat which differs from 
most other non-fibrous solids. The theoretical models based on this 
chain structure will be discussed later, as will the morphology of 
polyethylene. However, for now it should be noted that polyethylene 
is an example of a polymer which crystallizes partially, and thus is 
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neither completely crystalline nor completely amorphous. So, one 
parameter in any discussion of polyethylene must be of its crystal- 
linity, which is directly related to its density. Making the assump- 
tion that polyethylene consists of two distinct forms, the crystalline 
and the amorphous, it is much simpler, for this discussion, to des- 
cribe the specific heat of polyethylene in terms of the specific heat 
of these two idealized phases, i.e., the specific heat of any real 
sample can be expressed as an admixture of that of the two hypotheti- 
cal ones. In order to accomplish this, the heat capacity of various 
samples differing in crystallinity (density) was determined so that 
extrapolation to the completely amorphous and crystalline cases could 
be made. 

Heat capacity data for polyethylene is available from various 
sources in the temperature ranges from 1 to 5°K and from 20°K to its 
melting point (about 428°K) . The results of this research will 
bridge the existing temperature gap and should be of interest for two 
reasons. First, since polyethylene is an archtype of the chain 
polymers, the behavior of its heat capacity over the entire tempera- 
ture might be representative of polymers in general, and second, to 
allow a possible check on the aforementioned theoretical predictions. 

Previous Measurements 

Below 100°K there is not much experimental data on the specific 
heat of polyethylene. Above this temperature a considerable number 
of measurements have been made to determine the specific heat of the 
solid as it passes through the glass transition and the melting point. 
However these latter results are not germane to the present investiga- 
tion. 
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The data for the temperature range below 5°K were available from 
two sources. Isaacs and Garland reported measurements on a single 
high density sample in the 1.8-5. 3°K interval obtained by calorimetric 
methods [2l]. The heat capacity of three samples of various densities 
was measured by indirect means in the range of from 1 to 4.5°K by 
Reese and Tucker [34] . The density dependence noted before is quite 
apparent when all four sets of measurements are compared [$ 4 ] . It 
should be noted, however, that the data was analyzed assuming a cubic 
temperature dependence for the specific heat. This assumption is not 
verified by the present investigation, which shows a distinct non- 
cubic relation in this temperature region; however, the density de- 
pendence is valid. 

Dainton, Hoare, Evans, and Melia give results for two samples in 
the range 20-300°K [l3] . However no conclusions could be made about 

the density dependence, since the samples were very similar in density. 
The only other known source of specific heat data for polyethylene in 
this temperature range is the work of Sochava from 17 to 60°K [44]. 

The sample used by Sochava was not well characterized, the only infor- 
mation given being that it was a ’typical commercial product". Thus 
there is no way to determine the density relationships in the higher 
temperature range above 20°K. Above 90°K there are more measure- 
ments but it is felt that this information was not of interest insofar 
as the present investigation is concerned. 

Brief Description of Polyethylene 

Chemically, polyethylene is one of the simplest forms of chain 
polymers, consisting only of CH2 units. The bonding along the chain 
is covalent. Ihe chain length can be varied during manufacture, but 
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in general it is quite long, on the order of 1000 units. In addition, 
side branches can and do form on the chain, and is also a parameter 
which can be varied according to the manufacturing process. These 
branches, which are short chain sections normally only 2-5 units 
long, join to the main chain by the substitution of a carbon-carbon 
bond for a carbon-hydrogen bond [31,38] . 

Polyethylene is one of the polymers which crystallizes partially 
in the bulk form. A necessary condition for crystallization is geo- 
metrical regularity of the molecular structure, which is fulfilled by 
the simple form of the polyethylene chain [l0,3l]. This regularity 

hr 

allows the chains to form ordered arrays, which are considered the 
crystalline state. However on the basis of single crystals grown 
from solutions, which are quite small, it appears that the crystals 
are lamella, approximately 100 Angstrom units thick, with the molecu- 
lar chain axis being perpendicular to the face [23,35j . Whether this 
crystalline form is the same in the bulk is not known conclusively [ 31 ] 
However, it is well established that the crystalline part of poly- 
ethylene is composed of clusters of spherulites [ 31 J . 

The most important factor inhibiting crystal formation is branch- 
ing. This is evidenced by the fact that highly branched polyethylene 
displays relatively low crystalline content [.45] . This is quite 
reasonable since the regularity of the structure is broken at the 
branch points, thereby producing discontinuities or distortions in 
the crystal. 

If the chain is not in the well ordered array of the crystalline 
state, it is in some metastable position, intertwined and twisted 
around other chains. The complexity of the non- crystal line form makes 
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description very difficult. However, it sufficient for this discus- 
sion to characterize the non-crystalline part of the polyethylene as 
the amorphous regions. 



Since in the crystalline part of the polyethylene the molecules 
are more closely packed than in the amorphous, a higher density is 
found for the crystalline phase than for the amorphous [l5j . In 
fact if the density, ^ , of a particular sample is known, its weight 
percent crystallinity, X , can be calculated from the relation, 



1.1 






where and ^ refer to the densities of the crystalline and 

amorphous parts which are abtained by indirect methods [l5^ • The 
crystallinities of bulk polyethylene vary in general between 35 and 
95%. 

Heat Capacity Models 

The specific heat of a solid can be calculated if the tempera- 
ture dependence of its internal energy,^, is known. This follows from 
the definition 



1.2 
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If the internal energy’s temperature dependence is due to vibrations 
of the constituent parts of the solid about their equilibrium posi- 
tions, the total energy of oscillation over the whole range of the 
vibrational spectrum is, 



J f t)Jy~ , 
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where Zftyr)- 0 



is the average energy of a Planck 
oscillator with vibrational frequency v; and is the distribution 

function for the eigenf requencies of elastic vibrations of the solid, 
which is called the frequency spectrum ( 14] . If there are N mass 
points in the solid, then there can be only 3N normal modes of vibra- 
tion, assuming that the temperatures are sufficiently low that possible 
internal motion of the mass units can be neglected (the optical modes 
of vibration). This sets the restriction 

'6 



1.4 3/V - f . 



Rigorous treatment shows that the upper limits in the integrals in 
equations 1.3 and 1.4 can be replaced by a maximum or cut-off fre- 
quency [50]. Hence if jbr) is known, ^ can be determined from 
equation 1.4. 

This cut-off frequency can be interpreted as the frequency at 
which two neighboring mass units differ in phase by 180°. It should 
also be noted that is related to the forces acting between the 

neighbors, that is, greater forces imply larger cut-off frequencies, 
assuming particles of equal mass [ 14 ] . 

The central problem in the calculation of the internal energy 
of a system, and therefore its specific heat, is to find the frequency 
spectrum . This is in general extremely difficult. Therefore, 

several approximation methods have been employed for solids. One of 
the most successful approximations is the Debye model. 

The Debye approximations are that the frequency spectrum is that 
calculated assuming the solid is an isotropic continuum, subject only 
to the condition stated above concerning the cut-off frequency, and 
that there is a linear relation between the frequency and the wave 
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number, i.e., that there is no dispersion. These approximations can 
be expected to hold for long wavelengths (low temperature) [ 50 j . 
Debye theory than predicts a frequency distribution given by 

1.5 



o 






Using this in equations 1.2, 1.3 and 1.4 leads to the familiar Debye 
equation for the specific heat, 



1.6 



where 



CV « 3 R Pi f-r) 

P 3 ~e^i 



- 3 
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and ~&o ~ i- s a characteristic temperature of the material called 

/ 

the Debye temperature. For very low temperatures, 4.K ~^~o ’ 



equation 1.6 reduces to, 
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and for very high temperatures, 

1.8 C V = ZR T-*«- 

This last result is the same as is predicted classically using the 
equipartition of energy [ 24 ] . 

As a logical extention of this model, it is more reasonable to 
assume that there exists a cutoff wavelength, rather than a maximum 
frequency. Since it is known that for all solids the velocity of 
sound is dependent of the direction of propagation, the cut-off fre- 
quency will also be different for the different modes and directions 
[ 43 J. Recognizing this, equations 1.5 and 1.6 become, 
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where the subscripts refer to the transverse and longitudinal modes, 
and the )J^s are directly related to the average sound velocities in 
these two modes. It should be noted that this modification preserves 
the dependence of the specific heat at the lowest temperatures. 

Polymers differ from most other solids in the sense that they 
are essentially one dimensional in nature, because the forces acting 
along the chains are much greater than those between chains. If a 
Debye-type analysis is made for a one dimensional system, in a manner 
completely analogous to the preceeding, one obtains the following re- 
sults [ 47]; 
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The model employed above describes non-interacting chains, so 
is not of direct interest in this discussion. However, if the chains 
are intertwined and linked with one another, as in a chain polymer 
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such as polyethylene, there will be both the three dimensional and one 
dimensional continua existing together. This is the basis of a model 
due to Tarasov [ 47 ] . Since the forces acting along the chain are much 
greater than those between chains, it is expected that the maximum 
frequency for the one dimensional part of the model will be greater 
than that for the three dimensional part. Tarasov therefore predicted 
that there would exist two frequency spectra such that, 



where the subscripts refer to the one and three dimensional quantities. 
Using the Debye assumptions, frequency distributions can be found. 
However, the separation into transverse and longitudinal parts is 
again applicable, as it was above. Fbr the remainder of this deriva- 
tion only the longitudinal part will be considered, and the transverse 
contributions put in at the end by analogy with the Debye case. There- 
fore equation 1.13 becomes, 
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1.14 




for the longitudinal mode only (denoted by the superscript $ ). This 



equation can be rewritten as 
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where N 30 and N 13 refer to the first and second integrals above. 
The distribution of the modes are, for N ^3 
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If these spectra are the correct ones, the heat capacity is given by 
the resulting Tarasov equation, 
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for the longitudinal mode, where ~&i l } *& 3 l are related to 

and the other symbols have been previously defined. Inclusion of 

the transverse modes then yields the complete Tarasov relation, 

1.19 C v -Cy-t2C^. 



In the two temperature extremes, equation 1.19 reduces to 
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Between these extremes, there will be a region for which, 

+ *)] 
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No continuum model should be expected to explain the specific 
heat in the presence of dispersion, which will take place in any real 
lattice. In order to calculate the density of states in the vibra- 
tional spectrum correctly, the effects of all of the forces acting 
between the mass units must be included. In practice this is not 
feasible. As a result, the usual method of approach is to make real- 
istic estimates of the strength of the forces, and drop all but the 



dominate ones from the calculations. Blackman, and many others since, 



have achieved a certain amount of success along these lines; the de- 
gree of success depending on the validity of the assumptions about 
the relative strength of the forces, and the number of forces included 



ploys the discrete nature of the solid, should be more complete than 
the continuum treatment, if the correct parameters are used. There is 
however one feature of the lattice dynamics calculations which is the 
same as that obtained by the continuum assumptions: the heat capacity 
at the lowest temperatures should be proportional to the cube of the 
absolute temperature. 

For crystalline chain polymers an attempt has been made to obtain 
the vibrational spectrum using lattice dynamics by Stockmayer and 
Hecht whose work was greatly extended by Genensky and Newell fl9,46]. 
This analysis is based on the consideration that there are strong 
forces along the chains and weaker ones between them. The frequency 
spectrum deduced was quite complicated, but has certain characteris- 
tics which are predicted to be characteristic of the heat capacity of 
fibrous solids in the crystalline form, and thus should be applicable 
to crystalline polyethylene. These predictions include, (where 
now refers to the specific heat contribution from the interchain 
interaction, and C 3 to that along the chain) [ 19 J 



in the calculations [ 6,7,16,28]. This type of analysis, since it em- 
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where Ta* is proportional to the maximum vibrational frequency , 
oC is related to the nearest neighbor forces, and K to the forces which 
oppose the bending of the chain. The approximations are, 

1.29 P/*>2f , 0 /k > |0 , if- « 

where ^ is related to the force constant along the chain, and V is 
related to the second nearest neighbor forces. Using these relations 
temperature^ regions are predicted for which the heat capacity is at 
first (lowest temperatures) proportional to the cube of the tempera- 
ture, then changes to a T^ dependence, where b lies betweBn 2.5 and 
3, than b decreases to approximately .8 for T about one half T m . 

To complete this brief treatment of specific heat models which 
might be applicable to polyethylene, a related phenomenon, which 
appears in the low temperature heat capacity of some amorphous solids, 
must be mentioned. If the heat capacity of these substances is 
measured calor metrical ly at low temperatures, an "excess" specific 
heat is observed. To understand the term excess one should refer to 
equation 1.7 which can be rewritten as [24], 

r _ ^4! v — 

1.30 Cy g- r * 

where V is the volume of the sample, and ir is the velocity of sound 
at these temperatures, which can be interpeted as, 
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1.31 
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Therefore sound velocity measurements should yield the same value for 
C/T^ as those obtained from direct calorimetric methods. However, it 
has been found experimentally that for some amorphous substances, such 
as vitreous germanium, vitreous silica, and glycerol glass, the calor- 
imetric results are much greater than those given by the acoustic 
methods and are not T^ dependent [4,12,1?]. 

The explanation for this anomolous effect is that at low temper- 
atures there may be modes excited which are not measured by the 
acoustic methods. Since the mechanical sound measurements are con- 
cerned with the averaged (in the sense that the measurements must be 
macroscopic) elastic forces between the particles making up the solid, 
if there are low frequency modes which do not enter into these elastic 
vibrations, then the calorimetric measurements should be higher than 
predicted by equation 1.30. To explain the existence of these extra 
modes, Rosenstock has suggested that in disordered lattices, such as 
the amorphous solids, microscopic holes exist which may contain a 
particle which is bound in only one or two directions [39,40]. These 
loosely bound mass units can be likened to particles free to vibrate 
within a definite volume (the hole), and could thereby account for the 
low frequency non-acoustic modes. However, there are other hypotheses 
which attempt to explain the modes as being due to other phenomena 
such as non-symmetr ic lattice sites or distorted bonds [3,27]. In 
any case, the origin of these low frequency modes is not clear at 
this time. In the case of the examples cited above it has been possi- 
ble to empirically describe the low temperature departure of the 
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specific heat from the relation by adding a term (or terms) pro- 
portional to the Einstein specific heat to equation 1.30 [2,4,17] . 
This implicitly implys that some of the vibrating units do not con- 
tribute to the acoustic Sttodes and are accounted for in the total 
specific heat by an expression such as, 
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Yf is the excitation frequency for these modes, and Mf is the num- 
ber of vibrators involved [24]. Adding a term such as this to equation 
1.7 or 1.30 results in 

1.33 c v * A r 3 ■+ 8 t 



where A and 3 are constants. Since the specific heat is no longer 

o 

proportional to the cube of the temperature, a plot of the C/T versus 
temperature will demonstrate the presence of the excess heat capacity. 
2. EXPERIMENTAL DETAILS 

Introduction 

In all calorimetric investigations it is necessary to know two 
quantities; the temperature of the sample and the acount of heat added 
to it. In this experiment, meeting the former requirement was the 
more difficult. 

The temperature determination must be made quickly and accurately. 
Since resistance thermometry fulfills these requirements, doped ger- 
manium resistors were selected as thermometers. These were chosen 
rather than carbon resistors since germanium resistors are reputed 
to maintain their resistance-temperature relationship even after a 
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large number of cyclings between room temperatures and liquid helium 
temperature [l8,25] . Other experimenters have found that this is not 
true for carbon resistors, so that they must be recalibrated after 
each temperature cycling [l8] . The need for recalibration would 
greatly complicate the experiment thus making the choice of germanium 
resistors an easy one. 

Since the temperature range covered by this investigation was 
from about 2 to 30°K, doped germanium resistors (Minneapol is-Honey- 
well, models 2401 and 3401) were selected which had sufficiently high 
sensitivity, and yet, for the above temperatures, had resistance 
values within the range of ordinary laboratory resistance measuring 
devices. Once the resistors were chosen the most important task was 
to find the relation between the resistance and the absolute tempera- 
ture. The discussion of the calibration procedures is the subject of 
Appendix A. 

Having insured that a quick and accurate determination of the 
sample temperature and its changes could be made, attention was di- 
rected to the second requirement. Adding heat electrically to a 
resistance element in thermal contact with the sample provides a 
simple and accurate method of deliberately raising the sample's tem- 
perature. However great care must be taken to insurfi that this 
energy goes solely to the sample, and that it is the only source of 
heat input. To this end, an adiabatic calorimeter was constructed in 
which the sample could be suspended in vaccuo by high thermal resis- 
tance mountings from a heat shield, which completely surrounded the 
sample and was held at the same temperature. The design and construc- 
tion of this calorimeter will be discussed in a subsequent section. 
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Sample Preparation 



The heat capacity of three polyethylene samples of various den- 
sities were to be measured. Two of the samples were Marlex, manu- 
factured by Phillips Petrol ium Company, and were received in the form 
of extruded cylindrical pellets approximately 3 mm. in diameter and 
length. One of the samples, Marlex II, had branching introduced to 
impede crystallization. The third sample, a typical commercial low 
density polyethylene, was a 1" round bar. The samples are referred 
to as Marlex I, Marlex II, and Low Density in the remainder of this 
thesis. 

An important consideration in the preparation of the samples was 
that the time interval between the application of heat and the achieve- 
ment of thermal equilibrium be as short as possible so as to minimize 
errors due to heat leaks into or out of the sample. The low thermal 
conductivity of polyethylene at low temperatures. makes this a diffi- 
cult problem [34]. Because of the calorimeter design, the samples 
were to be cast into 9/16" rods about 6 cm. in length. If such a 
polyethylene rod were heated from the center, it was estimated that 
it would require times of about 100 seconds to achieve thermal equili- 
brium at 5°K and nearly 1 1/2 hours at 30°K. Such comparatively long 
thermal relaxation times would undoubtedly introduce serious errors 
in the determination of the amount of heat added to the sample unless 
exceptional precautions in the control of the temperature of the adia- 
batic shield were taken, in addition to making the experiment extreme- 
ly tedious to perform. Thus considerable effort was expended to pre- 
pare samples with shorter thermal equilibrium times. The method 
chosen was in inbed within the sample 99.999% pure copper wires such 
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that the wires were no farther than .3 cm. apart. In addition the 
sample was wrapped in a copper foil. The heat could thus be distri- 
buted throughout the sample by means of low thermal resistance copper 
paths. This was expected to greatly reduce the thermal relaxation 
times estimated above. The success of this method is seen by the fact 
that the time required during the actual experiment to achieve thermal 
equilibrium was never more than about 12 minutes at the highest temp- 
eratures . 

Introduction of the wires into the sample in a uniform fashion 
presented a considerable problem. The following system was eventually 
evolved. The pellets (or in the case of the low density sample, the 
shavings obtained from turning the bar on a lathe) were first cast 
into small rods of 1/4" and 3/16", and wrapped with 30" of the copper 
wire. The rods were compactly assembled within a 9/16" mold, and 
additional wires inserted into the resulting gaps. The sample was 
then cast at 155°C under 900 PSI for 16 hours, and then allowed to 
cool to room temperature (about 12 hours). 

The extra material added to the samples was to be considered 
part of the calorimeter’s sample holder. It consisted of 30 f of 
.0051" diameter drawn copper wire, .001" thick copper foil 5.580 by 
5.035 cm., and .08 grams of G.E. 7031 varnish. 

The polyethylene density was determined by measuring the density 
of the entire sample and correcting for the copper wire which was 
introduced. The crystallinity of each sample was then calculated by 
use of equation 1.1 with the values 
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where the and refer to pure amorphous and crystalline poly- 
ethylene densities ^9,36]. It must be noted here that the densities 
given above are not the only ones available. For the amorphous case 
there is fair agreement in the literature that .85 is a good average, 
but more attention must be given to the crystalline density [l5]. For 
this case, the density has been measured by X-ray methods to be as 
high as 1.014 gm/cm^ and as low as .964 gm/cm^ depending on the amount 
of branching present [ 49 j. Furthermore voids or free volume in the 
spherulitic crystalline form of bulk polyethylene could reduce the 
higher figure noted above even in the absence of branching [ 10,11 ]• 
These considerations will be referred to again when the error of the 
experiment is discussed. The measured densities and calculated crys- 
tallinities are given in the table below, where the error indicated 
in the latter quantity is that due to the uncertainty in the density 
measurement only. 

TABLE I 

Densities and crystallinities of the samples. 





Dens ity 


% Crystallinity 


Mar lex I 


.973 1 .002 gm/cm^ 


84.1 t 1.5% 


Marlex II 


.958 t .002 


74.5 1 1.6 


Low Density 


.915 1 .002 


45.9 t 2.5 



Calorimeter Design and Construction 



In order to accurately measure heat capacities, the sample must 
be thermally isolated to insure that the noted temperature differences 
of the sample are due entirely to the measured heat input. To accom- 
plish this, it is necessary to remove all conduction and radiation 
paths to and from the sample. The method employed here was to enclose 
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the sample in a shield which is kept at the same temperature as the 
sample, but insulated from it. 

In this experiment the heat capacity would be measured at various 
temperatures above the temperature of a liquid helium bath by success- 
ive applications of heat to the sample. Hence the shield could not be 
in good thermal contact with the bath. This thermal isolation was 
accomplished by suspending the shield inside of a vacuum container 
which was in turn in contact with the bath. A can type adiabatic 
shield was made from a copper 2" tube, 3” long, threaded at the top 
to mate with its lid. The lid was soldered to a 3/8" stainless tube 
(.01" wall thickness), which connected to a suitable vacuum container. 
The tube then extended out of the cryostat to peripheral vacuum equip- 
ment and served as the pumping line and electrical lead housing. The 
shield also had provisions for receiving a resistor which would sense 
the shield’s temperature. 

There would be a heat leak from the shield to the vacuum con- 
tainer through the support, but this was compensated for by means of 
an electrical heater on the shield. In addition this heater was used 
to raise the shield’s temperature to that of the sample. Seventeen 
feet of .002" manganin wire, coiled and glued with G.E. 7031 varnish 
to the top of the shield, served as the shield heater. 

The sample holder, in addition to holding the 9/16" sample, was 
used to mount the sample heater, temperature sensing resistor, and a 
receptacle which would be the contact point for the thermal switch 
(to be discussed later). A copper ring with the necessary provisions 
was constructed for this purpose. When heat capacity measurements 
were made, the measured values would be the total of the sample plus 
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the addenda (the holder, heater, thermometer, etc.). Although the 
latter’s heat capacity would be determined experimentally and sub- 
tracted from the total, it was still desirable to keep it as small 
as possible. Based on the size of the calorimeter, it was determined 
that a sample of 10 cm^ could be accommodated. With a sample of this 
size it was estimated (correctly, as was determined later) that the 
heat capacity of the addenda would amount to less than 10% of the 
total, which was deemed acceptable. 

The sample heater, which consisted of 15” of .002” manganin wire 
(100 ohms), was wou&d around the thermometer receptacle. A four wire 
connection was made to the heater, and the wires were thermally 
grounded to the adiabatic shield. The power was measured potentio- 
metrically using a Leeds and Northrup K-3 potentiometer. A Varian 
model 501A time mark generator, used in conjunction with a Hewlett- 
Packard model 521C electronic counter, determined the heating times 
(the timer was later changed to an Anadex model CF-200R counter- 
timer) and was triggered by the sample heater switch. 

It would be impossible, regardless of the care taken, to insure 
that the shield would be at exactly the same temperature as the sample 
at all times. However, if the shield were thermally isolated from 
the sample, small temperature differences could be tolerated. Based 
on this the sample holder was suspended within the shield by 6 nylon 
threads, approximately .12 mm. in diameter, and a minimum of 1.9 cm. 
in length. Assuming a maximum of .5° temperature difference between 
the shield and the sample, the heat leak through these supports would 
be on the order of 10“^ mj/sec which, for equilibrium times of less 
than an hour, could be neglected (radiation losses are of the same 
order of magnitude). 
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When the calorimeter was initially cooled, or when cooling back 
to the bath temperature after a series of measurements, a method had 
to be incorporated in the design to cool the isolated sample in a 
reasonable length of time. For this reason a mechanical switching 
device was constructed. This was chosen over the use of exchange gas 
in the vacuum spaces, since the latter method has been shown to have 
quite long recovery times and has attendent gas absorption problems 
[ 29 J. The switch consisted of a brass finger which would make contact 
with a receptacle on the sample holder. The finger was soldered to a 
brass rod which extended out of the vacuum container into the liquid 
helium bath through a nickel bellows and from there out of the cryo- 
stat, The finger could be raised and lowered manually so as to make 
thermal and mechanical contact with the holder whenever it was re- 
quired. The original design was modified after an initial trial by 
tinning the brass finger and its mating receptacle with metallic in- 
dium, for better thermal contact. In addition a brass leaf was added 
which would make contact with the adiabatic shield when the thermal 
switch was closed in order to cool the entire calorimeter simultane- 
ously. 

The two thermometers (resistors) were attached to the shield and 
the sample. Three wire connections, which were thermally anchored to 
the shield, were made to each resistor. The resistance measuring 
equipment was identical to that used by Reese and Tucker and will not 
be discussed here [ 34 J. 

The completed calorimeter is shown in figure 1. 
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Figure 1. Calorimeter with sample in place. 

a. Snmple 

b. Sample holder 

c. Sample thermometer 

d. Sample heater 

e . Thermal switch' 

f. Shield thermometer 

g. Nylon support 

h. Anchor 






Measurement Procedures 



In preparation for the heat capacity measurements, the sample 
was placed in the calorimeter and cooled to the desired initial tem- 
perature with the thermal switch closed. When the thermometers indi- 
cated that the sample and shield were at the liquid helium bath temp- 
erature, the thermal switch was then opened. This set up a small 
amount of frictional heating as the contact was broken, but within 
30 minutes the system would again be in equilibrium. 

The sample was then heated by applying electrical power to the 
sample heater in amounts which would raise the temperature of the 
sample in 5% increments of its average temperature. The power was 
recorded at each temperature step as well as the initial and final 
readings of the time and the thermometers. The heating times were 
kept between 20 and 60 seconds, which required that the power be in- 
creased every four to six measurements. When the desired final temp- 
erature was reached, the thermal switch was again closed, the sample 
cooled, and the procedure repeated. 

While the heat was being supplied to the sample, the adiabatic 
shield was also heated by means of its heater. This required care 

since, as previously mentioned, it was highly desirable that the 

✓ 

shield and the sample be at the same temperature. After several 
trials, however, certain regularities in the amount of heat needed 
for the shield at various temperatures became evident, and the desired 
result was obtained. It is felt that there never existed more than a 
one percent temperature difference between the two thermometers while 
actual measurements were being taken, and for the bulk of the runs, 
the difference was considerably less. 
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Due to the close proximity of sample thermometer to the heater, 
while the power was being supplied the temperature would be observed 
to overshoot the desired value, and, when the power was removed, drop 
back exponentially. The time required for the sample to reach equil- 
ibrium was from 2 seconds to about 12 minutes as the temperature was 
increased from its lowest value to its highest. These times were not 
considered excessive, except at the highest temperatures when the 
shield might tend to drift while equilibrium was being established. 
This drifting was countered by monitoring both thermometers and keep- 
ing the shield at the desired temperature by removing or increasing 
the shield power. 

The heat capacity of the samples was determined by five separate 
sets of measurements. The first three covered the range of 4.5 to 
30°K, one series for each of the three samples. The other two were 
for the temperature interval 2.4 to about 8°K for the Marlex I and 
Low Density samples only. For each set of measurements the aforemen- 
tioned temperature ranges were covered three times. 

Since the measured heat capacities were the sum of those of the 
sample and the addenda, the latter was measured in two separate ex- 
periments, for the temperature ranges of 4.5 to 30°K and 2.4 to 8°K, 
so that its contribution could be subtracted from the total. In each 
of the experiments, the addenda heat capacity was measured four times 
over the covered range. The procedures were identical to those des- 
cribed above except due to the reduced mass, the amount of heat re- 
quired and the equilibrium times were both substantially decreased. 

3. RESULTS OF THE EXPERIMENT 

Data Reduction 

When the experimental measurements of the heat capacity were 
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made, the data consisted of the power reading, the elapsed time inter 
val during which the power was applied, and the initial and final 
temperatures as indicated by resistance values- The heat capacity 
can then be calculated from, 



Due to the complexity of the resistance-temperature relation (Appen- 
dix A), the calculations were performed by the GDC 1604 digital com- 
puter. 

All of the data were reduced to heat capacities, and of the ones 
which displayed large deviations, only those attributable to obvious 
blunders in data recording were rejected. 

Since the heat capacities measurements for each sample represen- 
ted the total of the sample plus the addenda, the latter T s contribu- 
tion had to be taken into account and corrected for. As mentioned 
before, measurements of the heat capacity of the addenda were per- 
formed. These heat capacities were for discrete, random temperatures 
and, as such, were not useful for the adjustment which was to be made 
Therefore the data were best fitted, in the least squared sense, by 
the computer to a continuous heat capacity-temperature relation of 
the form, 



With this relationship the total heat capacity for each sample 
could be corrected so that the specific heat of the sample alone 
could be determined. These calculations were also performed by the 
computer and the resulting heat capacity data were again fitted to 



3.1 C 
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equation 3.2. The specific heat was obtained by dividing the heat 
capacity by the mass of the sample. 

The measured specific heat for each sample is found in Appendix 
C, and the parameters of equation 3.2 for the samples and the addenda 
in Appendix E. The smoothed results from the fitted curve are shown 
in Appendix D, and a specific heat-temperature graph is displayed in 
figure 2 where the difference between the smoothed curve and the actual 
data points is too small to be shown. 

The temperature range, total number of data points, and the per 
cent RMS deviation of the actual data from the smooth curve are pre- 
sented in Table II. 

Table II 

Temperature range, amount of data, and RMS fit of the data 
to a smoothed curve for each of the samples and the addenda. 





Temperature 

range 


Number of 
data points 


%RMS 

deviation 


Addenda 


2.4 - 


30° K 


218 


1.70% 


Marlex I 


2.4 - 


30 


161 


1.52 


Mar lex II 


4.5 - 


30 


117 


1.01 


Low Density 


2.4 - 


30 


158 


1.09 



Extrapolation to 100% and 0% Crystallinity 
As mentioned in the first chapter, the specific heat of any 
polyethylene sample might best be expressed as admixture of the specif- 
ic heats of the completely amorphous and completely crystalline phases. 
As suggested by the earlier low precision work of Reese and Tucker, ahd 
very well substantiated by the present investigation, the relation is 
linear [34], That is 

3.3 C = X C 10 o ■+ O'*) Co 
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Figure 2. Specific heat versus temperature. 
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where the C's on the right side refer to hypothetical values for the 
specific heat of crystalline and amorphous polyethylene. This lin- 
earity is reasonable if one considers that it is possible (in theory) 
to break a sample up into many very small macroscopic volumes, each 
of which is either entirely amorphous or crystalline (if the boundar- 
ies are neglected). The two idealized polyethylenes should have dif- 
ferent values for the specific heat at any particular low temperature** 
since the amorphous state should be characterized by a set of para- 
meters which will differ from those of the crystalline state inasmuch 
as the interchain binding scheme is completely different. 

Using this assumption, the heat capacities for the three sam- 
ples, divided by the cube of the temperature, were plotted as a 
function of crystallinity for various temperatures (Figure 3). The 
linearity is quite apparent, and on the basis of this, it is felt that 
equation 3.3 is indeed valid. Using this linear relationship one can 
perform a simple extrapolation to obtain values for the specific heat 
as a function of temperature for the hypothetical 100% and 0% cases. 
However, rather than use the C/T 3 relation, which is shown for des- 
criptive purposes only, the computer was used to perform the extrapo- 
lation from the equation 

3.4 C ~ w ^ 

where /vy\ and Jlr were determined from the experimental (smoothed) 
values of C, the specific heat for the three samples of differing 
crystallinities. The extrapolation to both 0% and 100% was performed 
every .04° and the results were fitted to equation 3.2, the same as 
the real samples. 

Appendix D presents the extrapolated data, and figure 2 is a 
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% CRYSTALLINITY 

Figure 3, C/T 3 versus crystallinity. 
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plot of the specific heat versus temperature. In Appendix E are 
shown the parameters for equation 3.2 for all of the real and extra- 
polated samples. 

The 100% Extrapolated Data 

The heat capacity obtained from the extrapolation to 100% 
crystallinity had some interesting features. A plot of C/T^ versus 
temperature shows a relatively large temperature region (up to about 
9°K) where the specific heat has a cubic temperature dependence (fig- 
ure 4). This is in agreement with the expectations if there exist 
relatively large binding forces in the crystalline case, since on the 
basis of the Debye model larger forces imply larger $ which in turn 
imply an increase in the temperature range in which the cubic depen- 
dence holds (equation 1.7). 

For this case, an excess heat capacity is not indicated by the 

3 

C/T plot, nor is one expected since this material should be a con- 
tinuous polycrystal. 

When the heat capacity for this hypothetical sample is plotted 
against temperature on log-log paper, it is possible to obtain the 
parameters A and a of the expression, 

3.5 C = AT*. 

There appeared two regions in which they were constant, indicating 
that the specific heat is proportional to some power, a, of the temp- 
erature. The first region was the cubic dependence noted before. The 
other was from about 20 to 26°K where the exponent, a, was approximate- 
ly 2.37-2.40. In between the value decreased uniformly. This is in 
disagreement with the prediction of Genensky and Newell who expected 
it to be equal to some constant slightly greater than 2.5 for an 
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C/T (MJ/°K -GM) 




Figure 4. 

A plot of C/T 3 versus temperature for polyethylene, 
where the solid lines are the experimental and the 
extrapolated data, the dotted line below the 0 % curve 
is the data minus the Einstein term, and the other 
dotted lines are the predicted T = 0°K values for 
these samples. The dots represent the predicted 
specific heat values in all cases. 
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appreciable temperature region, for a crystalline polymer [l9]. 

This contradiction will be discussed later. 

The Tarasov model could be well fitted to the extrapolated data. 
The method used to determine the parameters was essentially the same 
as that suggested by Wunderlich, except that the equations used are 
for the extended model, equation 1.19 rather than the one suggested by 
Tarasov, which doesn’t separate the longitudinal and transverse con- 
tributions [52]. According to equation 1.22 there should exist a 
region in which the 'lie's t capacity becomes linear in temperature. This 
will occur at a temperature well above that at which most of the three 
dimensional modes are filled and above the temperature range of this 
investigation. Therefore the data of Dainton, et al, giving the 
specific heat of polyethylene from 20-300°K, were used to determine 
the linear region [ 1 3 ] . This procedure will be justified later. The 
linear region appears to occur between 80 and 110°K, where C/T has 
the value 6 . 69m j/°K^-gm. Rewriting equation 1.22, it is seen that in 
this region, 



3.6 



c /t 




Another relation between and '&\ l . is necessary to determine them 

uniquely. This was accomplished by considering that the since it 
was the longitudinal mode limit for the chain, might be well character- 
ized by the largest frequency for which the C-C bond will be stretched. 
Reference to the Raman spectrum of polyethylene indicates a value of 
1150 cm*"l for this parameter, and therefore a value of 1655°K was 
assumed for 'O’il [ 15 ]. This fixed at 707°K. (The use of an appar- 
ent optical mode in this discussion, which up to now had only consid- 
ered acoustical vibrations, will be discussed later). 
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Having determined these two parameters, it is then necessary to 
fix the remaining two thetas needed for the extended Tarasov model, 
and^t . One relation comes from equation 1.20, which is here 
rewritten, 

3 7 C/7-3 =3 R ["T ^ 

This is the low temperature limit, where C/T^ is constant, and the re- 
sults below ^ 9°K could be used to evaluate the left hand side. In 
order to solve for either or another relation is needed. 

This was obtained from sound velocity measurements at liquid nitro- 
gen temperatures for a highly crystalline sample. These measurements 
were performed by Professor W. Reese [ 33 j. Extrapolation of the re- 
sults to zero degrees Kelvin predicted that the ratio of the longi- 
tudinal sound velocity to that of the transverse was 1.67. Since the 
thetas have the same relative relationship to one another as do the 
sound velocities, ^3r and -& 31 . could now be determined. The Tarasov 
parameters for this case are found in Table III. 

The 0% Extrapolated Data 

Extrapolation to the 0% crystallinity (completely amorphous) 
case led to results which were quite different from those discussed 
above. In particular, the plot of C/T 3 versus the temperature shows 
that the heat capacity was not proportional to the cube of the temp- 
erature anywhere. The curve starts with a positive slope at 2.5°K, 
increases to a maximum at around 5°K, then decreases uniformly. This 
is also the case for the real samples as can be seen in figure 4. 

Since even the Marlex I sample (crystallinity of 84.1%) displays this 
hump, while the 100 % extrapolated sample does not, it is quite apparent 
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that the presence of the hump must be due to the amorphous poly- 
ethylene. The hump is interpeted as an indication of an excess heat 
capacity. Referring to the first section the heat capacity for the 
lowest temperatures should then be of the form, 



C = AT 1 -V V £ ( J t') 

y= ZJkNe 



where A and are constants to be determined. This equation 
was best fitted, in the least squared sense, to the extrapolated re- 
sults (and to the Low Density data where a correction was applied to 
account for the crystalline part) to yield values for /, A and . 
For the model which describes the portion of the heat capacity not 
given by the Einstein term, one is forced to accept a continuum type, 
since it is not possible to do lattice dynamics for the amorphous 
case. From the foregoing discussion, the extended Tarasov model was 
chosen. Assuming this, A is interpeted as 

3.10 A-- - v) L “3 ( 4 I T + * 



There is no reason to assume that the one dimension parameters 
will change from those used previously, since they are merely func- 
tions of the forces along the chain. However another relation be- 
tween and is still needed in order to determine them uniquely. 

The ratio of the longitudinal and transverse sound velocities used 
earlier is not really applicable in this case, since the sound vel- 
ocities were those for a highly crystalline sample; but, lacking any 
other evidence, the same relation was assumed. 

Using the above assumptions, and the values of A and Y, the 
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Tarasov parameters can be found. They are presented in Table III 
along with oth&r results for this extrapolated hypothetical sample. 

Summary 

The percentage difference between the 0% and the 100% extrapola- 
ted samples were plotted against temperature as shown in figure 5. 
Although the curve is approaching zero, it is impossible to determine 
with any degree of certainty when if ever this would occur. However, 
it is felt that the percentage difference would be insignificant by 
100°K, where the data is linear in temperature [l3j. This supports 
the assumption that ‘d*, will be the same for all samples. 

The heat capacities for the real samples and the extrapolated 
ones are shown in Appendix D at .5°K increments, together with the 
specific heats predicted by the Tarasov equation, equation 1.19, which 
is here expanded as 




The superscripts on the refer to either the 0% or the 100% para- 
meter^, which are presented in Table III. The crystallinities, / , 
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are shown in Table I. 



Two hitherto undefined terms also appear in Table III. These 
are the Tarasov parameters which would result if no separation into 
the longitudinal and transverse modes were assumed, which is the ori- 
ginal model discussed by Tarasov. In this case, equation 3.11 would 
become, 



3.12 



C =3R{xfD,(%>^[D,(^)-D 3 (^)]] 






+ (i-y)XE(4)- 



In addition, the classic Debye temperature is shown as calculated 
from the low temperature limit. 

In figure 4 the dotted line beneath the 0% crystallinity sample 
curve is the result of subtracting the term from the data, to 

illustrate how it accounts for the excess heat capacity. For the 
other samples, the dotted lines indicate the 0°K intercept, and for 
all cases, the circled points are the values predicted from equation 
3.11. 

In Appendix D the smoothed values of the specific heat for each 
sample is shown, along with the calculated specific heat of the Tara- 
sov model as given by equation 3.11. This calculation was performed 
by the computer (the mass of the polyethylene vibrating unit is taken 
as 14.1 atomic mass units). 
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Table III 



The Debye, Tarasov, and excess heat capacity parameters, 
and the low temperature values of the specific heat for 
the extrapolated data. 







100% 


0% 




(mj/°K-gm) 


0.0 


3.24 


-3s 


<°K) 


260. 1 


175.4 


■Or, 


(°K) 


873.9 


873.9 




(°K) 


141.9 


78.6 




(°K) 


1655.0 


1655.0 


-&T 


(°K) 


707.0 


707.0 


G 3L 


(°K) 


222.8 


123.4 


3» T 


(°K) 


133.7 


74.1 




(°K) 


- 


23.0 


c -*£(30 

T3 


(mj/°K-gm, low 
temp. ) 


.00787 


.0256 



Accuracy of the Results 

To estimate the accuracy of the heat capacity data, four sources 
of error were considered; the determination of the amount of heat in- 
troduced into the sample, the temperature determination, the effect 
of the shield, and the heat capacity of the addenda. 

The amount of heat energy introduced into the sample was obtained 
from three measurements; voltage, current, and time. The voltages (and 
current) could be determined to six places on the potentiometer. Since 
the same power was used on from four to six successive measurements, 
it was possible to check the reproducibility of the readings. Only 
at the very highest temperatures, where the resistance of the heater 
was probably changing with temperature, did it vary as much as . 1%. 

The current measurement had the same inconsistency, but since it was 
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determined by measuring the voltage across a standard resistor, it 
also had an inherent maximum error of .05%. Timing was presumed ac- 
curate to one part in thirty thousand. The maximum total error from 
this source is therefore estimated to be .05% at the low end of the 
temperature sdale and .25% at the high end. 

As discussed in Appendix A, the thermometers were assumed cor- 
rect to within .5% of the absolute temperature. This uncertainty will 
influence the accuracy of the heat capacity results in two ways. Since 
a temperature difference is divided into the total heat to obtain the 
heat capacity, an error can arise there. However any experimental or 
systematic errors in the temperature determination will be the same 
over a short temperature interval and should cancel out (this is en- 
sured by the use of the smooth fitted curve which should average out 
these errors). Therefore the possible error from this source should 
be only that due to the error of judgment in reading the resistance 
bridge. This results in uncertainties in the heat capacity determi- 
nation of .12% at 2.4°K up to about .2% at 30°K. 

The other manner in which the inaccuracy of the temperature 
scale enters into the error is in assigning a temperature to a par- 
ticular heat capacity measurement. Since the specific heat can be 
expressed at any point as being proportional to some power of the 
temperature (the exponent will be between 3 and 2.4 for this case), 
the percentage error in the specific heat will be larger than that of 
the temperature. In this case the maximum possible error amounts to 
1.5% at the lowest temperature and 1.2% at the highest from this 
source. 

Since it was impossible to completely isolate the sample from 



48 



the shield thermally, considerable care had to be taken to ensure 
that they were at the same temperature, to prevent possible heat 
leakage. Experimentally, the effects of large temperature differ- 
ences were investigated, by purposely overshooting the shield tem- 
perature and observing the effect it had on heat capacity data. On 
the basis of these experiments it was decided that, with the experi- 
mental procedures used when actual data were being taken, this con- 
tribution to the error could be neglected for all but the highest 
temperatures, when, due to the insensitivity of the shield thermom- 
eter, difficulty insuring equal temperature of the shield and sample 
might occur. At these temperatures, the error might approach .5%. 

The heat capacity of the addenda was measured over the entire 
temperature range, so that it could be applied as a correction to 
the total measured for any sample. Since the mass of the addenda 
was chosen to be as small as possible, so that its heat capacity 
would be on the order of 10% of that of the samples, any possible 
errors which are present in the measurement of its heat capacity 
would be reduced by a factor of ten when applied to the total possi- 
ble error for a sample. Considering the preceeding this uncertainty 
amounts to approximately . 2% over the entire temperature range. 

Considering all of the above, it is felt that the measured heat 
capacities were accurate to within 1.9% for the temperature up to 
about 26°K, and 3% above that. 

The subject of Appendix B is the measurement of the specific 
heat of copper. The results are compared with data found in the 
literature which presumably is of higher precision and accuracy than 
this work. The comparison shows that the results are consistant, 
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within the error estimate, but the present results tend to be low 
over most of the temperature range, which could be an indication of 
a systematic error. This point was not investigated further. 

As mentioned before, there was no statistical criteria estab- 
lished for rejecting any of the data, and the only points discarded 
were ones in which a large deviation could be attributed to a blunder 
in data accrual. 

When the data were fitted to equation 3.2, the RMS deviation 
(Table II) was found to be well within the estimated accuracy of the 
experiment and therefore the use of the smoothed data is considered 
to be a valid representation of the results. 

Since one of the objectives of this investigation was to verify 
the density (crystallinity) dependence of the specific heat of poly- 
ethylene, and an extrapolation was made to the hypothetical 100% and 
0% crystalline phases, an estimate was made as to the reliability of 
this. The heat capacity of the three real samples would be exactly 
linear in crystallinity over the entire temperature range, if a varia- 
tion of at most 4% were made in any one of the heat capacities, or if 
a 1.5% adjustment were made in one of the crystallinities. Addition- 
ally there are other combinations of adjustments which could be made, 
but it is the uncertainty in the degree of crystallinity which is be- 
lieved to be subject to the most error. As noted in Table I, there 
is an inherent, approximate 2% error which is due solely to the un- 
certainties in the densities. In addition, since the densities of 
the amorphous and crystalline phases are not well established, it is 
estimated that the crystallinity of a sample could be in error as 
much as 4%, abased on the density method of determining crystallinities 
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(equation 1.1). In view of this, it is felt that the linearity is 
valid to the precision of this experiment. 

Due to the possible error in the crystallinities, the accuracy 
of the extrapolated samples can be better than 2% for the amorphous 
case, and . 87o for the crystalline phase. The uncertainties which are 
present in the actual specific heat is taken over into the extrapola- 
ted data. If the relative error in the heat capacity at any particular 
temperature is the same for all of the samples, as should be expected 
since they should exhibit the same systematic error if it is present, 
the maximum possible error in the extrapolated data due to this cause 
will be the same, or about 1.97o over most of the temperature range. 

Thus the uncertainty in these figures will be approximately 3 and 4% 
for the completely crystalline and amorphous phases. 

The precision of this experiment can be estimated by the close- 
ness of the fit of the actual. data points to the smoothed curve given 
by equation 3.2. .-.The deviations appear random and the per cent RMS 
deviation is given in Table II, where the % RMS deviation is defined 
by, 



3*13 



7oRMS = 100 % 




This precision, which is also an estimate of the reproducibility 
of the results, is between 1.01 and -1.52% for the polyethylene sam- 
ples . 

4. CONCLUSIONS 

Comparison with Other Results 

The results of this research have been compared with all of the 
known experimental work on the specific heat measurements on poly- 
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ethylene in their common temperature ranges. A detailed survey of 
these comparisons follows : 

(a) Isaacs and Garland [ 21 ] 

The specific heat measurements of Isaacs and Garland were on 
a single high density sample in the temperature range 1.8 to 5.3°K. 

They assumed that the specific heat in this range displaced T 3 depen- 
dence and therefore did not separately determine the heat capacity of 
the addenda, but rather calculated it and subtracted it from the total. 
This led as they report to an unaccountable deivation from the T be- 
havior, but they attributed it to experimental error. However it is 
felt that they were in reality observing the excess heat capacity 
which was disguised by their method of analysis. 

For the comparison in the mutual temperature range, it was noted 
that both the amorphous and crystalline parts of polyethylene are pro- 
portional to the cube of the temperature if the Einstein term is sub- 
tracted. Since their sample did not have the same degree of crystallin- 
ity of any of the ones measured in this work, the following equation 
was used to determine what specific heat should be expected on the 
basis of the present results; 




This is merely the low temperature limit of equation 3.12 with the 
measured values inserted. This assumed no model, however; only that 
the specific heat be T^ dependent at low temperatures, which it is 
for any realistic case, and that the linearity in crystallinity is 
valid. For the crystallinity of Isaacs and Garland’s sample, two 
values were investigated; first the reported X-ray determined result 
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of 90%, and second, a crystallinity of 85% deduced from the reported 
density. In Table IV where the comparison is made, the two different 
values for the C/T 3 resulting from the different crystallinities are 
subscripted to distinquish them. 

Examination of the table shows that the value deduced from this 
work with the density determined crystallinity is in excellent agree- 
ment with Isaacs and Garland’s, if it is recognized that they could 
not separate the excess heat capacity. Use of this crystallinity 
rather than the X-ray value contains a certain amount of consistency 
since the present results were based on the density measurements. It 
is noted they give no estimate of the possible error in the X-ray 
value, although it is doubted that it is as high as 6%. 

Table IV 



Comparison of present results with those of Isaacs and 
Garland- (Values expressed in mj/°K^-gm) 





C ^ r3 (I&G) 


C / T 3 (90%) 


C / T 3 (85%) 


T = 2°K 


.0113 


.0097 


.0108 


3 


,0113 


.0100 


.0112 


4 


.0113 


.0102 


.0115 


5 


.0113 


.0103 


.0116 



(b) Reese and Tucker [34] 

The results of thermal conductivity measurements on three sam- 
ples of polyethylene were used by Reese and Tucker to indirectly 
determine the specific heats of the samples. This led to a low pre- 
cision averaged set of values over the temperature range 1 to 4.5°K. 
The comparison with the present work was made in the same manner at 
it was for Isaacs and Garland’s work, and equation 4.1 was again used. 
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However in order to compare the values, a conversion was necessary. 

All of the results of Reese and Tucker were in terms of volume per- 
centages of crystallinities and the specific heats were in units of 
mj/°K 4 - cm3. The values used below are not as reported, but have been 
converted using their values of the densities. 

The present investigation was not independent of the work of 
Reese and Tucker, and, in fact, was an outgrowth of it. Therefore 
the same types of polyethylene were used in both experiments. How- 
ever slightly different densities (crystallinities) were found for 
the samples of the two investigations. This is explained by the dif- 
ferent molding processes, and cooling times, which is known to effect 
the degree of crystallinity [3l]. 

The comparison is shown in Table V. The sample names are from 
the work of Reese and Tucker and should not be confuted with the same 
names used in the present work. The calculated crystallinities are 
shown in the table. 

These results are in agreement within the reported error given 
by Reese and Tucker (10%) except for Marlex I. However in view of the 
indirect methods of their measurements, it is not considered exces- 
sive. In addition they report a "dish shape" characteristic if the 
values of C/T 3 is plotted against temperature. This was proposed to 
be due to a lack of precision in the temperature and thermal conductance 
determinations at the two temperature extremes of their range. If it 
is assumed that this was the case at only the low end, then the points 
above 2°K will display one half of the "dish" and will appear as the 
curves in figure 3 between 2 and 4.5°K. This is believed to have 
been the case since at the lowest temperatures, the parameter from 
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which the heat capacities was deduced (the thermal equilibrium times) 
would be greatly increased due to the thermal resistance at the sam- 
ple-low temperature sink interface. In light of the above, the 
results are considered to be in good agreement. 



Table V 

Comparison of present results with those of Reese and 
Tucker (Values expressed in mj/°K^-gm) 



T (°K) 

2 

3 

4 

5 



C/ r3 (R&T) 

Marl ex I (X = 7 9%) 
.0147 
.0147 
.0147 
.0147 



C/T^ (Present) 

.0118 

.0124 

.0128 

.0129 



2 

3 

4 

5 



Marl ex II (X = 62%) 



.0166 


.0149 


.0166 


.0159 


.0166 


.0167 


.0166 


.0168 



Low Density (X = 47%) 



2 

3 

4 

5 



0210 


.0175 


0210 


.0190 


0210 


.0200 


0210 


.0202 



(c) Sochava [ 44 ] 

As mentioned in the introduction, the specific heat measure- 
ments reported by Sochava for the temperature range 17 to 60°K are on 
a sample of unknown density. Since it has been shown here that the 
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specific heat below 30°K is strongly density dependent, comparison 
is difficult. However, if Sochava's data are placed on a specific 
heat-crystallinity plot from the present work and are adjusted to 
fall on the straight lines, a crystallinity of 53% is determined for 
Sochava's sample. This is a reasonable value for commercial poly- 
ethylene, and it is felt that, as far as the comparison can be made, 
there is agreement between the two sets of data. 

(d) Dainton, Evans, Hoare, and Melia [ 1 3 ] 

Dainton, et al , reported measurements on the specific heat 
of two samples of polyethylene for the temperature range 20 to 300°K. 
Their two samples were of very similar crystallinity (78 and 7 9%) , 
and due to the manner in which the data were presented, comparison 
could be made only at 20 and 30°K. Again the crystallinity was de- 
duced from the reported density (it was given for only one sample), 
but in this case it was identical with that given as determined by 
X-ray methods. To make the comparison predictions from the present 
results had to be adjusted to 78 or 79% crystallinities. This was 
accomplished by plotting the specific heat versus crystallinity for 
the three measured samples at 20 and 30°K, connecting the points with 
a straight line, and finding the specific heats corresponding bo the 
reported crystallinities. Table VT shows the results of this predic- 
tion, along with the results of Dainton 1 s et al . The subscripts 78% 
and 79% refer to the two samples. 
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Table VI 



Comparison of present results with those of Dainton, 



et al . 


(Values 


given in mj/°K-gm) 








Dainton, et al 


Present 


work 




c 78% 


c 79% 


c 78% 


C79% 


= 20°K 


66.8 


62.1 


62.2 


61.6 


30 


147.0 


146.2 


146.7 


146.4 



Considering that the 20°K results reported by Dainton, et al , are 
extrapolated from 22°K (the lowest temperature of actual measurements), 
the agreement is felt to be excellent. 

The Tarasov Model [ 47 ] 

The present results show clearly that a strong specific heat- 
crystallinity dependence is evident in polyethylene below 30°K. 

Hence no single set of Tarasov parameters can be expected to predict 
the heat capacity in general, and therefore a choice was made to sep- 
arate the contributions of the two phases, amorphous and crystalline, 
and obtain two sets of thetas. But they are not completely indepen- 
dent, as will be discussed. 

The density dependence does not extend much beyond the tempera- 
ture range of this investigation. This was determined from the fact 
that the results of Sochava and Dainton, et al , are essentially the 
same above 50°K, whereas they had differed by over 10% at 30°K. Addi- 
tionally data from Dainton, et al , on a sample of crystallinity 58% 
above 90°K show that the specific heat for this sample is the same as 
for the more crystalline ones. (Above 110°K this is no longer true, 
but it is due to effects which are of no concern here). This com- 
pares well with the indications of figure 5 which is a plot of the 
relative difference between the amorphous and crystalline results 
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Uo) 



(Cq - C 100 )/c 0 




Figure 5. Relative difference between the 0 $ and 
100$ extrapolated data versus temperature. 
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below 30°K. From the figure it is clear that the difference will be 
quite small for temperatures above about 60°K, if the curve is extra- 
polated above 30° K. 

The disappearance of the density dependence at the higher temper- 
atures has a rather natural explanation also on the basis of the 
Tarasov model . The three dimensional parameters associated with the 
two phases are small, which implies that the three dimensional modes 
will be nearly all filled at not very high temperatures. Thus the 
major contribution to the specific heat should be due to the one 
dimensional vibrations, which are assumed the same for both phases. 

The assumption that the one dimensional parameters are the same 
in both the amorphous and crystalline cases is based on the fact that 
these parameters are associated only with the forces along the chains 
which should be independent of the configuration the chain has assumed. 
It is also within the temperature range in which the density dependence 
disappears that the specific heat becomes linear in temperature, as is 
predicted by the Tarasov model, and it was from this temperature re- 
gion that the one dimensional thetas were determined. The linear 
portion of the specific heat-temperature curve yielded only an aver- 
age value of the thetas, and more information was needed to separate 
it into the longitudinal and transverse parts. 

The assumed value of 1655°K for the longitudinal one dimensional 
theta was based on the measured maximum frequency associated with 
the C-C stretch of the polyethylene chain as determined by Raman 
spectrum [ 15 ] . Since this parameter had been taken to be related to 
the maximum acoustic frequency along the chain, the use of what appears 
to be an optical mode may not be valid, but this seemed physically 
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representative of what this mode should be. In any event it should be 
pointed out that the predicted specific heat on the extended Tarasov 
model is relatively insensitive to variations of the values of 
and as long as the relation 



be .6 (the same ratio as assumed for the three dimensional case), it 
changes the specific heat preduct ion at 20°K by only .5%. 

One further remark must be made on the one dimensional para- 
meters. On the basis of experimental data, obtained by neutron 
diffraction, the upper limit of the bond stretching acoustic mode 
observed in polyethylene was found to be 545 t 25 cm“^ [ 32 J - This 
corresponds to a characteristic temperature of 785 t 35°K. If this 
value is assumed for the ^ , it predicts specific heats which are 
much too high from about 12°K on up in temperature. This is inter- 
peted not as anything wrong with the value of the maximum frequency 
noted above, but rather an indication that there is no real physical 
significance of the characteristic temperature from the continuum 



The three dimensional thetas for the amorphous and crystalline 
parts differ from one another in the manner in which they should be 
expected. They are related to the forces acting between the poly- 
ethylene chains, and, since the forces are most probably of the type 
which decrease rapidly with distance, the less dense and less closely 
packed chains of the amorphous state should have smaller thetas. The 
three dimensional parameter deduced from the experimental data support 



4.2 




is satisfied. For example, if the ratio of to is taken to 



model . 
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this qualitative description. 

The calculated values of the specific heat based on the original 
Tarasov model are not in good agreement with the experimental re- 
sults. As can be seen in Appendix D, there is a 5-6% difference 
around 20°K, although in other temperature ranges the fit is better. 
The extended Tarasov model, which separated the longitudinal and 
transverse contributions, was used since it gave an approximately 
2-4% improvement over the entire range when compared with the basic 
model. However in the final analysis the Tarasov model should not 
be expected to explain the specific heat of polyethylene any more 
than the Debye theory can for other solids, since they contain inher- 
ent errotteous assumptions concerning the lack of dispersion. In view 
of this, the fact that the Tarasov model fits as well as it does 
makes it useful as a first approximation for the specific heats of 
fibrous solids such as polyethylene at low temperatures. 

The Stockmayer-Hecht Model [ 19,46] 

In order to investigate the applicability of the model of Stock- 

mayer and Hecht as expanded by Genensky and Newell, the specific heat 

of the 100% crystallinity case was the only one considered, as the 

model was suggested for a fibrous crystalline solid. At the lowest 

temperatures, the model did predict the cubic temperature dependence, 

but the predicted temperature region for which the specific heat would 

2 5 + 

be proportional to T * was not found. However it is felt that this 
was an over simplification on the part of Genensky and Newell for, 
even if their analysis were correct, it would not be rigorously re- 
flected by this relation. The reason for this is that they predict 
in some temperature region the specific heat should be 
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C = At 



r/t 4 6T 3 

where A and B are constants (the first term is the contribution of 
the interchain interactions and the second is the result of vibrations 
along the chains). This prediction does not give any appreciable 
temperature region in which the specific heat is proportional to a 
constant power of the temperature. 

To test the model for the temperatures at which the specific 
heat was investigated, a plot of C/T^^ versus the square root of the 
temperature should have a region where the slope is constant, if the 
model is correct. Figure 6 shows that this does not occur in the 
temperature range below 30°K, and, since the slope is negative at 
the high temperature end, it is doubted that it will occur. 

Q / O 

However, as is also shown in figure 6, plbtting C/T versus T ' 
does show a linear region (22 to 27°K) which implys 

c/ n 

4.4 C - AT +ST . 

This could be interpeted as implying that their analysis is incorrect 

in the sense that they do not anticipate that the linear temperature 

contribution from the vibrations along the chain would occur at the 

5 /2 

same temperatures that the T relation holds for the interchain 
specific heat part. However, this cannot be reconciled with the model 
as presented, because the relative size of the two contributions is 
wrong. 

It is not possible in practice to draw from experimental heat 
capacity data reliable inferences regarding the detailed form of the 
frequency distribution, as the heat capacity tends to smooth out the 
individual features of the frequency distribution. Hence the compar- 
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C/t * A T 3 / 2 ♦ B 



30 60 90 120 150 

^3/ 2 ^ o 3/ 2 j 




Figure 6. C as two functions of the temperature 



ison with the model of Stockmayer and Hecht cannot specifically deter- 
mine where their analysis is in error. On the basis of the apparent 

relation for the interchain specific heat contribution (if the 
interpetat ion of figure 6 is correct), it seems that the essence of 
the model is correct, but that a reexamination of the one dimensional 
part is in order. However other assumptions of theirs may be wrong. 
Since the polyethylene crystal appears to be lamellae, (although not 
directly observed in the bulk polyethylene) and the crystalline form 
is spherulitic, the structure may be such that the idealized tetra- 
gonal crystal, assumed by Stockmayer and Hecht, will not yield the 
correct predictions. Moreover since they had to make estimates of 
the relative sizes of the force constants, of which little is known, 
erroneous assumptions here could lead to different forms for the fre- 
quency spectrum, and thereby errors in the deduced heat capacity. In 
any event in its present form the specific heat predicted from this 
model is not completely applicable to polyethylene below 30°K. 

The Excess Heat Capacity 

In the absence of any acoustical measurements made at low temp- 
eratures for the heat capacity of polyethylene, the humping of the 

o 

C/T versus temperature curve, figure 4, was interpeted as an excess 
heat capacity. The appearance of this excess was not expected, but 
it is not unreasonable that it does appear when it is recognized that 
it is present in other amorphous solids [4,12,1?]. The origin of the 
excess heat capacity i* unknown, but attributing it to the existence 
of non-acoust ical low frequency modes seems reasonable. In particu- 
lar, the Rosenstock theory is appealing since it is known that free 
volume (holes) does exist in amorphous polyethylene [lO,ll^]. However 
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no reasonable conclusions can be reached on this or any other model 
which might explain the excess on the basis of the present results. 

The addition of the Einstein term to the amorphous heat capacity 
empirically fits the data quite well at the low temperatures, and no 
other correction is deemed necessary to account for the excess in an 
analytic expression for the specific heat. In fact, it is felt that 
when low temperature acoustic measurements of the specific heat of 
polyethylene become available, the value for the extrapolated T = 0°K 
point will agree with the following relation which is based on the 
present results; 

4>5 c = (.n«- .on?x)r l . 

From the size of the constant which multiplies the Einstein 
function, it is readily deduced that of all of the possible vibration 
modes of the solid, only .17% are taken up in the excess term. This 
can be compared with the number of these modes in vitreous germanium 
(.4%) and vitreous silica (2%), but no conclusions can be drawn [4,17] 

Summary 

As has been shown, both the continuum, Tarasov, model and the 
lattice dynamical, Stockmayer-Hecht , model fail to completely predict 
the specific heat of polyethylene. This is not surprising as they 
are based on models which cannot be exact for differing reasons. The 
continuum model is at best a first approximation since it ignores the 
dispersion of the sound velocity which surely takes place in any real 
solid. Hence the spectrum of vibrational states is wrong. But, con- 
sidered as a simple two parameter approximation, the predictions are 
very good. The Stockmayer-Hecht model which attempts to account for 
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the dispersion fails for different reasons. Since this model is 
based on an assumption of the interparticle forces in various direc- 
tions, there are many parameters which determine the predictions. If 
some are chosen incorrectly, the whole model fails. On the basis of 
the present results no attempt was made to evaluate the parameters, 
since one of the basic predictions of the model did not materialize. 

The Einstein term adjustment for the excess heat capacity found 
in the amorphous polyethylene was adequate to explain the departure 
from the cubic temperature dependence of the specific heat at the 
lowest temperatures. But until such time that empirical evidence is 
brought forward to verify some model to explain it completely, it must 
remain as an a posteriori addition to explain the experimental re- 
sults. 

The hypothesis that the specific heat of polyethylene is linear 
in crystallinity was proven, it is felt, in the temperature range 
covered by this investigation, and extrapolation to the completely 
amorphous and completely crystalline cases is considered valid. It 
is also believed that the crystallinity dependence has disappeared by 
50-60°K (considering only the temperatures below 110°K) . 

The present low temperature specific heat measurements are more 
reliable than those given by the prior investigations (Isaacs and 
Garland, and Reese and Tucker), since the existence of the excess 
heat capacity is recognized. By extrapolating these results to 0°K 
and by using the method of interpolating between crystallinities 
noted before, this work can be mated to that of Dainton, et al, to 
give a continuous specific heat-temperature coverage up to liO°K, 
where many other data become available. The resulting information 
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should prove valuable to workers investigating models of crystalline 
and non-crystalline fibrous solids into which nearly all polymers can 
be classed. 
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APPENDIX A- Calibration of the Thermometers 



Introduction 

Probably the most exacting aspect of this experiment was the 
determination of the absolute temperature scale. The difficulty 
stems from the fact that between 5 and 14° K there is no convenient 
temperature standard. Below 5°K the scale had been set by inter- 
national convention in terms of the vapor pressure of liquid helium, 
and between 14 and 25°K, the vapor pressure of liquid hydrogen is a 
good standard [8,48,5l]. Additionally the internationally accepted 
platinum resistance thermometers can be used above 10°K [20]. Be- 
tween 5 and 14°K the temperature scale is normally established by 
means of a gas bulb thermometer. 

In this investigation the resistance thermometers were calibra- 
ted below 4.2°K by means of a vapor bulb thermometer, and from 4.2 
to 30°K by means of a gas bulb thermometer. Liquid hydrogen vapor 
pressure thermometery was not attempted due to its inherent dangers: 
the facilities were not available to obviate them. Initially it was 
planned that the platinum resistance scale was to be compared with 
the gas bulb results, but was not since another comparison became 
available (the specific heat of copper, Appendix B), and it soon ap- 
peared that the higher accuracy in temperature range above 10°K was 
not necessary for the purpose of this investigation. 

The theory, construction and experimental procedures concerned 
with the calibration of the resistance thermometers is the subject of 
the remainder of this appendix, together with an estimate of the 
error in the temperature determination. 
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Theory of Constant Volume Gas Bulb Thermometry 
The equation of state for an ideal gas is, 

a. l = kVL 

TV r t 

where a and / refer to the initial and final states of the system. 

If the volume is held constant and the initial conditions are known, 
then there exists a linear relation between the temperature of the 
gas and its pressure. A constant volume gas bulb thermometer could 
consist of a large volume, , at the temperature to be measured, 
and a pressure sensing device at temperature,""^ , with volume ^con- 
nected by a tube of small volume V* 1 . Then equation A.l must be modi- 
fied to account for the fact that not all the gas is at the low 
temperature in the final condition; some is still in the gauge and 
connecting tube. To account for this equation A^l is replaced by, 



A - 2 'rfa* *+*)*?(* + ■% 

where Py is the pressure of the system when it is at Ty , and P is 
the pressure when the bulb is at temperature T. The integral ex- 
presses the fact that the temperature along the connecting tube 
varies from T to 7y . However, it is convenient to assign a temper- 
ature Tc such that, 



t' r v 



A. 3 



T c = <r‘ f 

t*T 



&l(r *( t ) 

t 



Solving equation A. 2 for"T, it is found that, 



A. 4 



T - T v P 



V 



P v fr-f *-) - Ptr+tr'lf 9 ' y-PZL 



) 
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Consider now the case where ^^0, that is assuming that the 
connecting tube’s volume is quite small in comparison with that of 
the other components. Then equation A. 4 becomes, 



A. 5 



T = 



PyO+tf-P'? 

Then if / , and if the bulb is at low temperature, most of the 
gas will be in bulb, and the pressure will be essentially linear in 
temper at ure. 

However, in an actual experiment the gas will not be ideal, as 
implied in the previous arguments. The equation of state in this 
case can be written, 



A. 6 



Zsf- »/?(';+ 4 •") 



T ' v ’ ’ ^ 

where the B± f s are the virial coefficients, and IT is the specific 

volume. If helium were used in the temperature range of the present 

experiment 8 y is of 2 or 3 orders of magnitude smaller than 8/ , and 

the term containing B z can be neglected [22]. Since (3, is also small, 

/?T 

IS can be set equal to ~p~ in the term which corrects for the second 
virial coefficient, so that, 



A . 7 pv = A, /?(Tf BP) 

putting -8 . 

Using A. 7 equation A. 2 becomes, 

V 



A ' 8 f> (r^e v p, + T.-H3.P, ' + r v ^ v Py) P (t->bp 1 1 T If Bj 



V 4 ~' r ' 



where the subscripts on the refer to the temperature at which 
they are evaluated. 

Finally, consider the form that equation A. 8 would take if the 
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initial conditions were varied somewhat, that is, by filling the 
system with the gas when the temperature of the bulb is at a third 
temperature T* and therefore has volume ^ due to thermal contrac- 
tion. Under these circumstances, equation A. 8 will be, 

- +— £ — )=. 

T v *£ y /» v T c *8 e p„J 

-y— -+ — -f tr' - \ 

T-iBP ^ T y +B y P ^ T c '+t3 c 'P) 

where the Pv refers to the pressure in the system when the bulb is at 
temperature TZ • (Notice that the T c's on the two sides of the equa- 
tion are not the same since the lower temperature is different.) 

Construction Gas Bulb Thermometer 

As indicated in the comment following equat ion fii. 5 , temperature- 
pressure linearity can be obtained if the gauge volume to bulb volume 
ratio is made small. As a design criteria was taken as 1/20, and 

the connecting tube's volume (T* would be held small compared to the 
volume of the gauge. 

To obtain the desired temperature range, the bulb was initially 
cooled to 4.2°K by a liquid helium bath, and then heated electrically 
above this point. The bulb was made of copper, to insure uniformity 
of temperature throughout it, and the resistors were placed directly 
on the outer surface of the bulb. 

The pressure measuring device was a Wallace and Tierman gauge, 
type FA-145 (2 turn dial, 0-800 mm. Hg range). A .44 mm. inside 
diameter stainless steel capillary tube served as the connecting 
tube, and the gas used was helium. 

The gauge was calibrated against a mercury manometer (assumed 
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correct) both before and after calibration of the resistors. There 
appeared a linear deviation with pressure which was regarded as a 
systematic error in the gauge and a correction for it was made in the 
subsequent calculations. 

The volume of the gauge, was determined by comparison with a 
known volume. The value was 

L r s 10.4(1 + 1.4 X 10“ 4 P) Z . 3 cm 3 

where P is the pressure in mm. of Hg. (This correction term accounts 
for the fact that the gauge had a bellows type pressure sensing ele- 
ment) . 

Since it was desired that V the volume of the bulb should be 
at least 230 cm^. A 1 1/2 M copper pipe, 10 M long, sealed at the ends 
produced a volume calculated to be 27 3.4 1 .7 cm^ at room temperature. 

Because of the size of the dewar, it was necessary to have the 
capillary 48" long. This resulted in a calculated volume for the 
connecting tube of only .19 cm^, well within the design criteria. 

As liquid helium temperature (approximately 4.2°K) would be the 
lowest temperature used in the calibration, some means must be de- 
vised to thermally isolate the bulb from the bath, once it was in 
equilibrium with the bath, in order to raise it above this point. 

For this reason the bulb was suspended in a vacuum container which 
was to be immersed in the bath. The bulb support was a 3/8" O.D. 

(.01" wall thickness) stainless steel tube. This tube also extended 
up and out of the cryostat and served as the pumping line. An elect- 
rical heater of about 1500 ohms, consisting of 12 feet of .002" man- 
ganin wire ( ^ = 78.3 ohms/ft.) was wound around the bulb and glued 
there with G. E. 7031 varnish. The heater was necessary to raise the 
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temperature and to compensate for the resultant heat leak through the 
support (and to a much lesser extent, the heat leak due to radiation 
and conduction by any residual gas). 

The resistors to be calibrated were inserted into small recep- 
tacles in the bulb. In order to avoid the relatively large electrical 
lead resistance, a three wire connection was made to each resistor. 
This, and a discussion of the low power A. C. bridge used, is given in 
the work by Reese and Tucker, and will not be presented here [34]. 

A schematic diagram of the completed system in shown in figure 7. 

Gas Bulb Calibration Procedures 

The gas bulb was filled with gaseous helium to about 790 mm. of 
Hg when the system was at liquid air temperature. The vacuum spaces 
were then pump&d to approximately 10“^ Torr and the liquid helium 
space was filled. 

The data accrual was quite simple. The procedure was to heat 
for a period of about 30 seconds, then turn off the heating power, and 
then observe the resistances come to equilibrium, which indicated 
temperature equilibrium. The pressures and resistance values were 
then recorded. The tendency for the temperature to drift downward, 
due to the heat leak, was compensated for by the heater. 

The observations were taken over a pressure range of about 40- 
300 mm. of Hg, corresponding to temperatures of from 4 to 30 degrees. 
They were taken twice over the entire range, at about 2-4 mm. steps, 
to check reproduct ibility, and a third time over the lowest portion 
of the range where the resistors were the most sensitive. Seventy- 
three points were taken altogether. 

After completion of this procedure, the bulb was allowed to come 
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To manometer 




a. Vapor bulb 

b. Resistors to be calibrated 



Figure 7. Schematic of gas bulb thermometer. 
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to complete equilibrium with the liquid helium bath. Using the vapor 
pressure, the bath temperature was determined. This would later 
serve as a check on the accuracy of the calibrating procedures. The 
helium was then evaporated from the dewar and the space filled with 
liquid nitrogen (99.9% pure). The pressure in the gas bulb together 
with the known temperature of the nitrogen bath gave values of 
and X/ which could then be used in equation A. 9, and the pressures 
could be converted into temperatures. 

Equation A. 9 is here rewritten with a slight amplification of 
the terms. 

fry* , ir , P« cr ' 

A * 10 TstiS^p* Tv-tS y l°„ T u -tf3 c p„ 

PV Ptr— 

T-+ RP T v ifSyP ^ -rj + filP 

where is the pressure in the system when the bulb is at liquid 
nitrogen temperature, is the liquid nitrogen temperature, V is 

the volume of the bulb at temperatures less than 30°K, and V* is the 
volume of the bulb at Ty . The subscript Y refers to ambient condi- 
tions in the laboratory, and the other symbols have the same meaning 
as before. Using published values for the thermal contraction of cop- 
per at the various temperatures, it is found that V r 270.8 1 .7 cm^ and 
V v s 27 3.0 1 .7 cm^ £26]. For an estimate of the X S , a value of 145°K 
was assumed foT 1 c , and 108°K for . These were obtained by assuming 
a linear temperature distribution along the capillary. This is suf- 
ficiently accurate as this part of the system is quite small. Values 
of the second virial coefficients were taken from Keeson;£22] 
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A. 11 



0 = (.2453 - 6.2009/T) °K/mm of Hg 
= .000163 
Q y = .000162 

8 C = .000179 
&c =.000186 . 

Using the above in equation A. 10 leads to an equation quadratic in 
Solving that for "T yields, 

A. 12 T- ■+ xfjf- 




r-Pl — - l. <r' _ \ _ n ( Lr tr' 

N \T»+e„p>j r Tt+n^J r \~r v t3 v p~ h TjfBc'P 

fJL -- ^ * L.IOO^ , 

Equation A. 12 was programmed for the GDC 1604 digital computer. 

The procedure discussed above gave only a pressure-resistance 

relationship until equation A. 12 was solved. To obtain the desired 

temperature-resistance dependence, all the data was best fitted, in 

the least squares sense, by the computer to the following type of 

equation for each resistor; 

1a - a - J 

A. 13 T = £ A: (Ary,. R) 

/ * I 

where the A*' $ and the best fitted sr\ (integers between 1 and 7) 
were determined. The values are presented at the end of this appen- 
dix, and the RMS deviation from this fitting is calculated to be .23%. 
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Calibration of the Thermometers below 4.2°K 

For temperature determination from 1 to 5.2°K, the International 
Commission on Weights and Measures adapted the "1958 He^ Vapor Pres- 
sure Scale of Temperatures" as the international standard [8,48]. 
Tables based on this scale give a conversion from vapor pressures to 
temperatures, hence the temperature can be ascertained from a reading 
of the vapor pressure alone. 

For this lower temperature determination, it was decided to use 
a vapor bulb thermometer, rather than merely measuring the pressure 
above the bath since the latter procedure involves making corrections 
and assumptions of questionable validity. 

Vapor bulb thermometery is performed by placing a container into 
the helium bath at the liquid level at which the temperature is de- 
sired to be known, and filling the container with gaseous helium 
until the gas condenses which indicates equilibrium between the 
liquid and gas phases at this pressure [42]. The pressure is then 
observed, and tables give the temperature directly. 

For the actual construction of the bulb a small copper cylinder 

o 

with a volume of approximately 2 cm J was connected by a 1/8 f stainless 
tube to a mercury manometer and filling manifold. The bulb was 
mounted directly on top of the vacuum container, and connecting tube 
was brought outside of the cryostat enclosed in an evacuated tube (to 
prevent gas oscillations in the system [42]. See figure 7. 

For the calibration, exchange gas was introduced into the vacuum 
space where the resistors were located, and the pressure above the 
bath was lowered in steps. The manometer readings were taken and 
the corresponding resistances were recorded. In all 32 points were 
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taken in the temperature range of from 2.4 to 4.2°K. 

After conversion of the pressures into temperatures, the resis- 
tance values were related to the temperature by equation A. 13. 

Recalibration of the Thermometers 

Inasmuch as the comparison experiment (Appendix B) indicated that 
a systematic error might be present in the temperature scale, the 
scale was reestablished. Also, in the initial calibrations, the vapor 
pressure and the gas bulb temperature scales (below and above 4.2°K) 
were independently determined with the result that there existed a 
slight discrepancy at their common temperatures. (Approximately .2% 
difference). Since it was desired to have a continuous scale over the 
entire range, the calibration was repeated insuring that the tempera- 
ture overlap was consistant with both types of temperature determina- 
tion. This was accomplished by using the liquid helium temperatures, 
as determined by vapor pressures, as the fixed points for the gas bulb 
calibration rather than the liquid nitrogen point. As a check as to 
the precision of this method, the predicted liquid nitrogen tempera- 
ture was only . 3°K different from the assumed value of 77.35°K, a 
difference of only .23%. 

The equipment and procedures were the same as outlined earlier. 

A total of 99 points were taken over the entire range, varying the 

i 

intervals so that a uniform coverage was obtained. The gauge volume 
was again experimentally determined, and was within 3% of the original 
measurement. The gauge was likewise recalibrated against a mercury 
manometer twice. 

The resistance values were related to the temperature by use of 
equation A. 13. 
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Accuracy of the Calibration 

In the temperature range below 4.2°K where the vapor bulb ther- 
mometer was employed, the only significant source of error in the 
temperature determination was the reading of the manometer. It could 
be read to within .001 cm of Hg with reproducibility of .003 cm. of 
Hg. This results in a maximum uncertainty in the temperature for this 
range of 3 mill idegrees . 

For the gas bulb calibration, which covered the temperature range 
of from 4.2 to 30°K, the following effects were considered to ascer- 
tain an estimate of the error; 

(a) The fixed temperature point. 

The calibration over the entire range is very sensitive to 
any errors in the determination of the fixed temperature point which 
essentially sets a value for the right hand side of equation A. 12. To 
determine this parameter, six independent observations of the tempera- 
ture were made, by vapor bulb methods, at the liquid helium bath tem- 
perature, and the pressure associated with each recorded. The fixed 
point was taken to be the average of these data. The resulting error 
in the temperature determination over the entire range due to this 
was estimated to be at most .003 degrees. 

(b) Errors in the measurement of volumes. 

Since the bulb’s volume was much greater than the other vol- 
umes in the system, the temperature is essentially linear in the 
pressure. This tends to make any error in the measurements of the 
volumes self compensating, and their net contribution to the error in 
the calibration is negligible. 
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(c) Thermalmolecular correction. 

Following Roberts and Sydor iak, the thermalmolecular cor- 
rection is .004° at the low end of the temperature scale (4.2°K), 
and negligible above 7°K [37]. This systematic error could have been 
corrected for, but was not in view of the existing larger errors which 
will be discussed below. 

(d) Pressure measurement. 

The pressure measuring gauge readings were compared with those 
of a mercury manometer both before and after the resistor calibration. 
This comparison was made stepwise going up in pressures, since this 
was the way data were taken, to include any hysteresis of the gauge’s 
bellows as part of the correction. The gauge was determined to have 
a linear error as a function of pressure, that is, the deviation was 
consistently high at the higher pressures, and low at the lower ones. 
The data were corrected for the presence of this systematic error. 

The statistical uncertainty of this correction is entirely compatible 
with the precision to which the gauge could be read, about .1 mm. Due 
to the uncertainty in the pressure determination, a possible error of 
.01° exists throughout the temperature range. 

Other sources of error, such as those considered by Beattie, and 
those due to the uncertainty of the values of the virial coefficients, 
accuracy of the mercury manometer, etc., were considered to be insig- 
nificant compared with the ones discussed above [l5]. 

Once the temperature scale was established, there exists the 
necessity of relating the resistance values of the resistors with the 
temperatures. This introduces another possible reading error which 
amounts to .1% over the entire range. However it is felt that due to 



84 



the insensitivity of the resistors at the higher temperatures, a more 
realistic value of .3% should be assigned there. After the resistance 
temperature relation had been established by equation A. 13, it was pos 
sible to assign the maximum temperature differences these errors would 
introduce, .0015° at 2.5 and .03° at the highest temperatures. It is 
also noted that to use the resistors as thermometer, another reading 
error of this size could exist. 

Summarizing all of the above, and recognizing that the fitted 
relation between the temperature and the sample thermometer resis- 
tance values contains an uncertainty, it is felt that the temperature 
determination by use of the resistors was accurate to within .5% of 
the absolute temperature. 

As an estimate of the reproducibility and precision of the cali- 
bration, the RMS deviation of the actual points from the fitted 
resistance-temperature was used. This amounted to .23%, and appeared 
to be randomized. Hence the statistical scatter of the calibration 
points from the fitted relation was well within the estimated errors 
of the procedure. 
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Parameters for the temperature-resistance fitted relation 



7 - £ 

/<- / 

Thermometer 





Shield 
(4. 2-30°K) 


Shield 
(2. 4-4. 20 K) 


Sample 
(2. 4-30° K) 


n 


6 


6 


4 


Al 


. 81986X10 -1 


-.11574X10- 1 


-.24787X10* 


A2 


. 65400X10 -2 


. 66904X10 -1 


. 32888X10 


a 3 


50749X10 1 


. 25718X10 -1 


-. 10029X10' 


A 4 


. 17961X10 2 


.10939X10° 


-. 36394X10 : 


a 5 


. 14366X10 2 


- . 27080X10 1 


.25491X10 


A6 


15249X10 3 


-. 69149X10 1 


-.45361X10 


a 7 


•19880X10 3 


. 38145X10 2 


. 36772X10 


% RMS 

Deviation .20 


.00 


.23 



86 



APPENDIX B - Comparison Experiment 

There has been considerable discussion about the need for, and 
choice of, a suitable calorimetric standard in the temperature range 
below 20°K [53]. Because of the large number of specific heat deter- 
minations for copper, its available purity, and other reasons, there 
is much support for its choice as such a standard. At the same time 
that this experiment was being performed, Martin reported the specific 
heat of a cast 99.999% pure copper sample in the temperature range of 
this experiment [30 ]. He asserts that his results are accurate to 
better than 1%, which implys a temperature determination to within 
.35% accuracy. 

To ascertain an estimate on the accuracy of the temperature 
scale used in the present experiment, it was decided to duplicate 
Martin's measurements. A single crystal copper sample of 99.9995% 
purity was obtained (Semi-Elements, Inc.), and was etched with nitric 
acid to a size which was compatible with the calorimeter. The re- 
sulting sample was a 2%" by 9/16" rod with a mass of 83.181 grams. 

Heat capacity measurements were performed on this sample, in the same 
manner as for the polyethylene samples, over the temperature range 
4.5 to 30°K. 

In performing these measurements, the only experimental procedu- 
ral differences from the polyethylene runs were that the amount of 
heat introduced for the temperature steps was considerably less, and 
the equilibrium times were much shorter, due to the differing heat 
capacity and thermal conductance. Three complete sets of measure- 
ments were obtained over the aforementioned temperature range, re- 
sulting in a total of 119 data points. 
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After smoothing the data by least squares fitting it to a seven 
parameter equation (equation 3,2) and correcting for the heat capacity 
of the addenda, the results were compared with those of Martin. They 
matched fairly well at the low temperature end, this work being about 
2% lower, and excellent agreement was obtained at the highest tempera- 
tures where no difference was found. However Martin’s results were 
over 5% higher at the midpoint, 17°K. The differences along the scale 
were not random but showed a clear divergence up to a maximum, and then 
converged to zero at the highest temperatures. This type of deviation 
can most easily be explained as some sort of systematic error in the 
determination of either temperature scale, or possibly a drift of the 
resistance values of the thermometers. Martin discounts this latter 
possibility in his work by repeated checks against both vapor bulb and 
gas bulb thermometery . In the present work nearly every time the sys- 
tem was at equilibrium with the helium bath, the thermometers were 
checked against the vapor pressure temperature of the bath. On the 
basis of this, it is felt that there was no resistance drift of the 
thermometers used here at least at the low temperatures. 

The size of the maximum difference in the two sets of results 
is quite large, but it should be pointed out that at 17°K an error 
in the temperature determination of only 1.3% (.2°) would explain 
the difference. This value of error is outside of the estimated 
error in the temperature scale used here, however, and so it was de- 
cided that a recalibration of the thermometers would be performed at 
the conclusion of the specific heat measurements on polyethylene 
(the copper data were taken about half way through the polyethylene 
work). It was felt that with the low temperature techniques acquired 
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by this experimenter that a more precise calibration might be achieved. 
The recalibration is discussed in Appendix A. 

After the recalibration was performed, the heat capacity data 
for copper were recalculated, and again compared with the results of 
Martin, The same relative deviation persisted, however it was approx- 
imately halved at the 17°K point (2.6% low). At the highest temper- 
atures investigated, 30°K, the results were again the same as Martin’s, 
but the same 2% difference was still present at 5°K. 

The new calibration was decided to be the more accurate and was 
used based on the following considerations. First the agreement with 
Martin’s copper data was better. Next the recalibration was per- 
formed after this experimenter had considerably more experience with 
the low temperature techniques, and was able to look for small sys- 
tematic errors which may have been overlooked the first time. Lastly 
this calibration was continuous over the range 2.4 to 30°K which was 
highly desirable (see Appendix A). 

At the conclusion of this work, Alhers published specific heat 
data on copper over essentially the same temperature range as that 
under consideration [ ij. His work compares very well with that of 
Martin, but has the same relatively deviation from this work. For 
comparison, the present copper specific heat data are presented be- 
low along with that of both Martin and Alhers. Additionally the 
parameters for the fitted curve equation (3.2) are given, along with 
its o/o RMS deviation from the actual data points, (.69%). 
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Parameters for the specific heat-temperature fitted relation 



9 * / 

c = £ A-T a ~' (C in mj/°K-gm) . For T between 2.4 and 30°K. 

A- i 

A 1 -.12463X10° 

A 2 . 96209X10 -1 

A 3 23948X10 -1 

A 4 .41835X10-2 

A 5 -.26648X10-3., 

A 6 • .10 385X10" 4 

A? -.14265X10-6 

% RMS .69 

Deviation 

Molecular weight of copper 63.57 gms/mole 
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COMPARISON Ol- SPECIFIC HEAT OF COPPER 





SPECIFIC HEATS 


IN MJ/OK/GM 




TEMPER AT uRE 


C ( TUCKER ) 


CIMART IN) 


C( ALH6RS) 


5.0 


.1443 


.1475 


. 1484 


6 . 0 


.2229 


. 2265 


.2275 


7.0 


^3283 


.3334 


.3339 


8.0 


. 4657 


. 4729 


.4725 


9.0 


. 6403 


.6505 


.6488 


10.0 


.8573 


.8718 


.8690 - 


11.0 


1.1226 


1.1444 


1 .1399 


12 . 0 


1 .4430 


1.4737 


1 . 4690 


13.0 


1.8266 


1.8708 


1.8647 


14 . 0 


2.2823 


2.3409 


2.3356 


15 . 0 


2.82o5 


2 . 896q 


2.8910 


l6.n 


3.4521 


3 . 54 4 o 


3 . 54 o o 


1 7 . 0 


4.1892 


4 . 3 o 0 6 


4.2917 


18 . 0 


5.0439 


5.1619 


5.1548 


19 . 0 


6.o286 


6 . 1542 


6 . 1388 


2 0 • 0 


7.1546 


7.2829 


7.2549 


21 . 0 


8.4324 


8.5537 


8.5193 


22.0 


9 , 8 7 o 2 


9 . 9695 _ 


+ 


23.0 


11.4734 


11.5367 


♦ 


24.0 


13.2436 


13.2685 


+ 


25. Q__ 


15.1772 


15 . 1584 


♦ 


26.0 


17.2649 


17.1997 


* 


27.0 


19.4898 


19.3727 


♦ 


28.0 


21.8264 


21.6576 


+ 


29.0 


24.2389 


24 . 0413 


* 


30 . 0 


26.6797 


26.5238 


* 




* INDICATES 


NO DATA 
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APPENDIX C - Raw Data 



MARLEX 1 



SPECIFIC HEATS 


IN MJ/OK/GM 


TEMPERATURE 


C 


2.597 


.18796 


2.604 


.19134 


2.587 


.19166 


2.741 


.23133 


2.768 


. 23560 


2.754 


.23643 


2.8 96 


. 27448 


2.913 


.28187 


2.934 


.28523 


3.048 


.32351 


3.0T4 


. 33065 


3.098 


.33663 


3.200 


.37506 


3.234 


.38555 


3.258 


.39376 


3.379 


.44463 


— 3.433 


.45855 


3.415 


.46099 


3.594 


.52877 


• 3T6T8 


. 52958 


3.617 


.53548 


3.808 


.61506 


3.804 


.63022 


3.819 


.63708 


3.998 


.71218 


4TC0 3 


.72960 


4.046 


.75817 


4.196 


.82776 


4T210 


.85138 


4.271 


.88733 


4.326 


.92379 


4.340 


. 93608 


4.366 


.96433 


4.4 19 


.98215 


4.490 


1.04819 


4.517 


1 .05959 


4.539 


1.07225 


4 • 6 37 


1 .13612 


4.635 


1 .14261 


4.673 


1 .16034 


"4.696' 


1 .18922 


4.774 


1.21758 


4.824 


1.29738 


4.845 


1.30108 


4.863 


1.32905 
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SPECIFIC HEATS IN MJ/OK/GM 



TEMPERATURE C 



4.970 
-4x99 r 
5.011 
5.10G 
5.135 
5.210 
5.229 
-5.304 
5.355 
5.418 
"5.468 
5.520 
5.604 
-5x6-37 
5.686 
5.746 
-5x771- 
5.91 1 
5.930 
6 ^086 
5.992 
6.181 
6.233 
6.236 
6.513 
"6x5 32- 
6.549 
6.572 
6.868 
6.905 
6.953 
*>.892 
7.178 
7.204 
-7.229 
7.344 
7.456 
7.535 
7.541 
7.711 
-7.84 7- 
7.909 
7.944 
0x152 
8.208 



1 .39622 
1T44958 
1 .45996 
1.54250 
1.61394 
1.64436 
1.66338 
1.74434 
1.79458 
1.84677 
1.91434 
1 .96849 
2.06463 
-2.08688- 
2.16203 
2.21610 
2.27986- 
2.42215 
2.42923 
2.62558 
2.62725 
2.75895 
2.81262“ 
2.84967 
3.19788 
3.24379- 
3.30149 
3.31813 
-3.72186- 
3.82183 
3.85513 
3.89666 
4.22355 
4.23521 
4.32056- 
4.50184 
4.63985 
4.83053- 
4.85262 
4.97886 
-5.41547- 
5.51952 
5.53639 
-5x94872 — 
6.12060 
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MARLEX 1 



SPECIFIC HEATS 


IN MJ/OK/GM 


’ TEMPERATURE 


C 


8.277 


6.29253 


8.437 


6.53970 


8.671 


■ 7.00062 


8.623 


7.19002 


8*773 


7.29157 


8.985 


7.70739 


9.075 


7.98402 


9 . 07 7 


8.034 1 1 


9.353 


8.54090 


9.415 


8.88594 


9.424 


8.93725 


9.610 


9.34961 


9.820 


9.74374 


9.847 


9.97924 


r 9.990 


10.61 803 


10.265 


1 1.22138 


10*33 1 


1 1 . 48649 


10.604 


12.30049 


10.834 


12.99738 


TOTS 92 


13.1 7960 


11.148 


14.03310 


11.339 


14.67437 


IT. 4 00 


1 4.93354 


1 1 .644 


15.82448 


1 1 .786 


16.29977 


11.949 


16.91263 


12.162 


17.62685 


12.198 


1 7.79657 


12.832 


20.28792 


12.859 


20.48844 


13.239 


22.24001 


1 37642 


23.724S5 


13.635 


23.74350 


13.856 


24.79318 


1 4.326 


26.70627 


14.343 


26.93295 


14.563 


27.92274 


14.988 


29.89960 


14.973 


31 .75482 


15.371 


31.86838 


~ 1 5.634 


32.10680 ~ - 


16.320 


37.06861 


17.374 


42.64875 


1?Y488~ 


43.70005 — 


17.819 


45.38054 
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MARLEX 1 



SPECIFIC HEATS IN MJ/OK/GM 



TEMPER A TIRE" 



C~ 



13.430 
-13.4 66 
18.875 
19.274 
-19.555 
19.799 
20.105 
-20.724 r 
20.884 
21.133 
-21-.667 
21.910 
22.468 
-23.02a 
22.828 
23.499 
-23.8 45 
24.186 
24.980 
25.001 
25.292 
26.191 
26.578 
27.159 
27.572 
28v825 



48.82613 
-49.36476 
' 51.37767 
54.50185 
-56.01417 
57.71516 
59.66975 
63.97106 
65.13815 
68.59212 
71.51811 
71.80872 
77.27470 
80.64122 
81.65808 
85.04005 
88.50399 - 
91.55677 
95.85445 
-98. 14700 
100.16850 
106.23490 
115.46458 
118.91252 
121.90134 
132.66855 
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MARLEX 2 



SPECIFIC HEATS IN MJ/OK/GM 



TEMPERATURE 


C 


4.349 


1 .10836 


4 .3T5 - 


1.1 5563 


4.373 


' 1.15943 


4.538 


1.26999 


4.540 


1 .28241 


4.568 


1.31700 


4.703 


1.41884 


4.74u 


1 • 46093 


4.762 


1.50724 


4.881 


1.60343 


4.948“ 


1 .66051 


4.967 


1 . 70458 


5.092 


1.82800 


5.1T0 


1 .91897 


5.187 


1.96224 


5.301 


2.05350 


S • 3 9 0 


2 . 17302 


5.484 


2.22102 


5.412 


2.25201 


51 6TF6 


2.44090 


5.670 


2.51030 


5.645 


2.53578 


5.850 


2.77896 


5.859 


• 2.78156 


5.897 


2.89658 


6*078 


*“■ 3.09759 


6.197 


3.28717 


6*165 


3.31436 


6.378 


3.60245 


6.454 • 


3.75170 


6.595 


3.93908 


6.718 


4.14388 


6.770 


4.29967 


6.952 


4.57285 


7.030 


4.73857 


7.118 


4.93654 


7.315 


5.26753 


7.330 


5.29696 


7.489 


5.69568 


7.642 


5.90489 


7.698 


6.04339 


7.875 


6.53743 


7.990 


6.69330 


0.098 


' 6.94474 


8.338 


7.47455 
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MARLEX 2 



SPECIFIC HEATS IN MJ/OK/GM 



TEMPERATURE 



C 



8.612 
—8.-673- 
8.720 
9.1 64 
-9.-337 
9.447 
9.71 1 
-9.907- 
10.000 

10.314 
-10.456- 

10.534 
10.860 
-1 1 . 0 17 - 
11.098 
11.391 
-11 .600" 
1 1.694 
1 1 .938 
12.1 95- 

12.315 
12.493 
12.797 
12.961 
13.066 

-13.4 13- 
13.701 
13.660 
-V4-.-1 44- 
14.353 
14.422 
-14.929- 
15.030 
15.236 
-15.738 
15.784 
16.060 
- 16 . 573 - 
16.596 
16.911 
-17.397 
17.448 
17.827 
- 1 8 . 1 45 
18.345 



8.36381 
8.37944 
8.42132 
9.63058 
10.13203 
10.62664 
1 1.61380 

-1 1 .85851 

12.28690 

13.08954 

-13.67324 

14.10426 

15.03178 

-15.61562 

16.19879 

16.97032 

17.80159 

18.33617 

19.13534 

20.23688 

20.92201 

21.41557 

22.72404 

23.71420 

23.89595 

25.43271 

26.84314 

26.95160 

29.02634 

30.01182 

30.75287 

-32.89604- 

33.22668 

34.73176 

-37.14855— 

37.24379 

39.29161 

-41.48706 

42.02919 

44.06952 

46.24928 

46.91189 

49.84552 

50.75218 

52.51535 
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MARLEX 2 



SPECIFIC HEATS IN MJ/OK/GM 



TEMPERATURE 


C 


18.832 


55.59494 


1 8.789 


5d .66262 


19.311 


' 58.92028 


19.688 


62.09741 


1 9 .790 


62.71279 


20.324 


66.22665 


20.554 


68.59343 


2 0 • G 3 3 


71 . 12623 


21.507 


73.83693 


21.511 


75.11742 


— — 21 .9 38 


78.73929 


22.729 


82.89937 


22.647 


84.30054 


23.086 


88.27922 


23.871 


92.26345 


23.844 


96.20608 


24.266 


96.79725 


24.965 


102.41275 


25.217 


104.56977 


25.469 


1 10.322^0 


26.1 18 


1 12.90491 


26.716 


117.49972 


26.707 


119.71567 


27.999 


131.33693 


28.399 


131.63365 


28. 6 or 


1 32 • 40289 


29.244 


137.68534 
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LOW DENSITY 



SPECIFIC HEATS IN MJ/OK/GM 
TEMPER A TUVE C - 



2.521 

■2.535 

2.579 

2.621 

2-.65S 

2.695 

2.731 

■2T78CT 

2.816 

2.851 
■2790a 
2.9U 1 
2.978 
-3.06-5 
3.097 
3.088 
-3.241 
3.251 
3.266 
'3'. 4 5 4" 
3.449 
3.440 
-3.645- 
3.648 
3.657 
-3.84-3 

3.851 
3.863 
4v052 
4.061 
4.064 
4v2 7 8 
4.328 
4.318 
4v335' 
4.474 
4.498 
4.497 
4.508 
4.642 
4". 682 
4.706 
4.831 
4v888 
4.932 



.28781 

“729686 

.31867 

.33245 

.35445 

.37404 

.38825 

“741757 

.42907 

.44835 

—.4771 1 

.50078 
.51675 
-.571 80- 
.57700 
.59378 
.65976— 
.67721 
.69384 
—81479 
.81595 
.82114 
-.96564 — 
.96623 
.97461 
1 .-13362- 
1.14367 
1 .14469 
1 .33429 
1.33820 
1 .34516 
1 .57952 
1.61775 
1.62234 
1.63967 
1.81021 
1.81792 
1.84558 — ' 
1 .86155 
2.03621 
-2-.06809- 
2.13112 
2.30063 

^.38673 

2.46734 
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LOW DENSITY 



t 



SPECIFIC HEATS 


IN MJ/OK/GM 




T ffMPER AT OR E 


C 




5.034 


2.63293 




5.116 


2.74904 




5.140 


' 2.79630 




5.189 


2.88276 




5.245 " 


2 . 96659 




5.360 


3.19051 




5.382 


3.21864 




5.450 


3.36224 




5.470 


3.36955 




5.615 


3.64384 




5.634 


3 . 68241 




5.712 


3.83U01 




5.713 


3.85035 




s.B9i 


4.22231 




5.916 


4.27282 




5.980 


4.32430 




5.992 


4 . 45900 




6.209 


4.91476 




6.235 


4.95730 




6*2 f 7 


5.03374 




6.303 


5.10324 




6.550 


5.69680 




6.575 


5.76155 




6.592 


5.78959 




6.644 


6.00095 




6T904 


6.59770 




6.928 


6.63617 




6.922 


6.65371 




7*268 


7.52172 




7.263 


7.56141 




7.287 


7.61844 




7.625 


_ 8.57451 




7.669 


8.71615 




8.024 


9.80153 




8.082 


9.98707 




8.405 


1 1.06141 




8.460 


1 1.21426 




8.505 


11.36980 




8.815 


12.46119 




8.899 


12.75976 




8.933 


12.86687 




9.231 


1 3.96903 




9.361 


14.44571 




9 . 3S7 


14.58159 




9.661 


15.61894 
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LOW DEMSITY 



SPECIFIC HEATS IN MJ/OK/GM 

temp era Tore C 



9.884 
9.894" 
10.127 
10.464 
10.452 
10.659 
11.038 
TT-.064- 
11.194 
1 1 .631 
TTt6'6'9 
1 1 .788 
12.229 
T 2.276 
12.410 
12.845 
T 2.9 01- 
13.063 
13.515 
T3t576- 
13.767 
14.245 
14.31 1 
14.520 
15.022 
T 5.094 
15.340 
15.862 
-15.940- 
16.220 
16.750 
T 6 . 8 4 5 
17.137 
17.700 

17.785 
18.713 
18.777 
19.002 

19.786 
19.860 
-19.986- 
20.904 
20.958 
2T.054- 
22.047 



16.52142 
-16.5793T 
1 7.59983 
18.91642 
10.95116 
19.73098 
21 .5751 1 
21 .69968 
22.29614 
24.39704 
24.42918 
25.18032 
27.30692 
27.65680" 
28.33392 
30.68352 
50.86788 
51.55914 
34.28889 
34.61570 
35.76025 
38.38456 
38.66656 
39.96585 
42.96486 
43.32534 
44.96633 
48.09940 
48.4751 1 
50.30252 
53.82994 
54.36106 
56.28382 
59.96677 
60.70527 
67.33329 
67.69524 
69.52506 
75.32322 
76.15117 
76.79421- 
84.22005 
84.50975 
85.05182 
92.92935 



/ 
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LOW DENSITY 



SPECIFIC HEATS 


IN MJ/OK/GM 


~TEMPERATCRF“ 


C 


22.1 1 8 


94.18537 


22.11 8 


94.42999 


23.228 


103.33531 


23.336 


103.43852 


^3*21 1 


103.96273 


24.385 


113.53377 


24.352 


115.23389 


24.560 


116. 08967 


25.490 


124.87124 


25.557 


125.00688 


25.776 


1 25.96857 


26.735 


133.83623 


26.654 


134.54346 


26.992 


135.1 4483 


27.869 


140.81310 


27.939 


143.93052 


28.21 5 


145.08745 


29.122 


154.61557 


29.283 


156.00964 


29.038 


161*87336 



i 



i 

r 

I 

1 
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THE 100 0/0 CRYSTALINE SAMPLE 



HEAT 


CAPACITIES 


IN MJ/OK/GM 




TEMPERATURE 


C ACTUAL 


TARASOV C 


PER CENT 








DIFFERENCE 


2.5 


.124 


.123 


.9 


3.0 


.211 


.213 


-.7 


3.5 


.335 


.337 


-.9 


4.0 


.501 


. 504 


-.6 


4.5 


.715 


.717 


-.3 


5.0 


.984 


.984 


.0 


5. 5 


1.313 


1.310 


.3 


6.0 


1 .708 


1.700 


.5 


6.5 


2.174 


2.162 


.5 


7.0 


2.715 


2.700 


.6 


7. 5 


3.337 


3.321 


. 5 


8.0 


4.04 5 


4.050 


.4 


8.5 


4.84 2 


4.832 


.2 


9.0 


5.732 


5.735 


-. 1 


9. 5 


6.719 


6.742 


-.3 


10.0 


7.805 


7.858 


-.7 


10.5 


8.994 


9.087 


-1.0 


11.0 


10.288 


10.433 


-1 .4 


11.5 


1 1.689 


1 1 .898 


-1.8 


12.0 


13.197 


> 13.485 


-2.2 


12. 5 


14.815 


15.195 


-2.6 


13.0 


16.543 


17.029 


-2.9 


13. 5 


18.381 


18.985 


— 3 . 3 


14.0 


20.330 


21.064 


— 3.6 


14. 5 


22.589 


23.263 


-3.9 


15.0 


' 24.558 


25.581 


-4.2 


15. 5 


26.836 


26.013 


-4.4 


16.0 


29.221 


30.557 


-4.6 


16.5 


31.713 


33.210 


-4.7 


17.0 


34.310 


35.966 


-4.8 


17.5 


37.010 


38.322 


-4.9 


18.0 


39.812 


41.773 


-4.9 


18.5 


42.714 


44.815 


-4.9 


19.0 


45.715 


47.942 


-4.9 


19.5 


43.813 


51.151 


-4.8 


20.0 


52.007 


54.435 


-4.7 


20.5 


55.295 


57.792 


-4.5 


21.0 


58.677 


61.215 


-4.3 


21.5 


62.152 


64.701 


-4.1 


22.0 


65.720 


68.245 


-3.8 


22.5 


69.381 


71.844 


-3.5 


23.0 


73.136 


75.492 


-3.2 


23. 5 


76.990 


79.187 


-2.9 


24.0 


80.942 


82.924 


-2.4 


24. 5 


84.996 ' 


86.701 


-2.0 


25.0 


89.156 


90.513 


-1.5 


25.5 


93.429 


94.358 


-1.0 


26.0 


97.821 


98.233 


-.4 


26.5 


102.339 


102.135 


.2 


27.0 


106.994 


106.061 


.9 


50.0 


332.800 


290.367 


12.8 


100.0 


682.400 


644.550 


5.5 


150. 0 


942.900 


920.149 


2.4 


200.0 


1197.000 


1 117.31 3 


6.7 
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THE MARLEX ONE SAMPLE 



HEAT CAPACITIES IN MJ/OK/GM 



TEMPERATURE 


C ACTUAL 


TARASOV C 


PER CENT 
DIFFERENCE 


2.5 


.172 


. 1 71 


.5 


5 . 0 


.30 3 


.303 


-.0 


5 . 5 


.487 , 


.490 


- .5 


4.0 


.755 


.739 


-.5 


4.5 


1 . 05 ? 


1 .057 


-.3 


5.0 


1.449 


1.448 


.1 


5 . 5 


1.920 


1.920 


.4 


6.0 


2 .4.95 


2.478 


.7 


6.5 


• 3.154 


3.127 


.8 


7.0 


3.909 


3.875 


.9 


7.5 


4.761 


4.725 


.0 


8.0 


5.714 


5.682 


.6 


8.5 


6.768 


6.749 


.3 


9.0 


7.926 


7.932 


-.1 


9 . 5 


9 . 187 


9.231 


-.5 


1 G . 0 


1 0 .553 


10.647 


-.9 


10.5 


12.023 


12.184 


- 1.3 


n .o 


13.597 


13.841 


- 1.3 


11 . 5 


15.275 


15.618 


- 2.2 


12.0 


17.057 


17.515 


- 2.7 


12.5 


18.943 


19.530 


- 3 . 1 


13.0 


20.931 


21.663 


- 3.5 


13 . 5 


23.021 


23.911 


- 3.9 


14.0 


25.213 


26.272 


- 4.2 


14.5 


27.506 


28.742 


- 4.5 


15.0 


29.899 


31.319 


- 4.7 


15.3 


32.392 


33.998 


- 5.0 


16.0 


34.985 ■ 


36.776 


- 5.1 


16.5 


37.676 


39.648 


• c. 


17.0 


40.465 


42.611 


- 5.3 


17.5 


43.353 


45.659 


- 5.3 


18.0 


46.338 


48.789 


- 5.3 


18.5 


49. 420 


51.995 


- 5.2 


19.0 


52.600 


55.274 


-5 . i 


19 . 5 


55.876 


55.621 


- 4.9 


20.0 


59.248 


62.031 


- 4.7 


20.5 


62.717 


65.501 


-4.4 


21.0 


66.281 


69.026 


-4.1 


21 . 5 


69.940 


72.602 


-3.8 


22.0 


73.693 


76.226 


- 3.4 


22.5 


77.539 


79.893 


- 3.0 


23.0 


81.476 


83.601 


- 2.6 


23 . 5 


85.505 


87.346 


- 2.2 


24.0 


89.622 


91.125 


- 1.7 


24 . 5 


93.825 


94.934 


-1 .2 


25.0 


98.112 


98.772 


-.7 


25 . 5 


102.479 


102.634 


-.2 


26.0 


106.923 


1 06.520 


.4 


26.5 


1 1 1.439 


110.426 


.9 


27.0 


1 16.021 


114.350 


1 .4 


50.0 


332.800 


296.206 


1 1 .0 


100.0 


682.400 


646.873 


5.2 


150.0 


942.900 


921.362 


2.3 


200.0 


1197.000 


1 118.064 


6.6 
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THE MARLEX TWO SAMPLE 



HEAT CAPACITIES IN MJ/OK/GM 



TEMPERATURE 


C ACTUAL 


TARASOV C 


PER CENT 
DIFFERENCE 


2.5 


* 


.201 


* 


3.0 


« 


.357 


« 


3. 5 


a 


.582 


* 


4.0 


* 


.831 


* 


4.5 


1 .249 


1.261 


-1.0 


5.0 


1.735 


1.729 


.4 


5.5 


2.315 


2.288 


1.2 


6.0 


2.995 


2.947 


1.6 


6.5 


3.778 


3.711 


1.8 


7.0 


4.666 


4.584 


1.8 


7. 5 


5.662 


5.573 


1.6 


8.0 


6.765 


6.679 


1.3 


8.5 


7.977 


7.907 


.9 


9.0 


9.297 


9.258 


.4 


9. 5 


10.724 


10.733 


-.1 


10.0 


12.258 


12.352 


-.6 


10. 5 


13.898 


14.054 


-1 . 1 


11.0 


15.643 


15.398 


-1 .6 


11.5 


17. 49'^ 


17.863 


-2.1 


12.0 


19.440 


19.947 


-2.6 


12.5 


21.491 


22.148 


-3. 1 


13.0 


23.642 


24.462 


-3.5 


1 3. 5 


25.891 


26.886 


-3.8 


14.0 


28.240 


29.417 


-4.2 


14.5 


30.685 


32.051 


-4.4 


15.0 


33.229 


34.784 


-4.7 


1 5. 5 


35.869 


37.612 


-4.9 


16.0 


38.606 


40.531 


-5.0 


16.5 


41.439 


43.536 


-5.1 


17.0 


44.370 


46.623 


-5.1 


1 7. 5 


47.398 


49.787 


— 5.0 


18.0 


50.522 


53.024 


-5.0 


18.5 


53.744 


56.330 


-4.0 


19.0 


57.062 


59.701 


-4.6 


19.5 


60 .476 


63.131 


-4.4 


20.0 


63.985 


66.617 


-4.1 


20.5 


67.588 


70.155 


-3.8 


21.0 


71.283 


73.742 


-3.4 


21.5 


75.067 


77. 372 


-3.1 


22.0 


78.956 


81 .044 


-2.7 


22.5 


82.887 


84.753 


-2.3 


23.0 


86.913 


88.497 


-1.8 


23. 5 


91.007 


92.272 


-1.4 


24.0 


95.161 


96.076 


-1.0 


24.5 


99.365 


99.905 


-.5 


25.0 


105.607 


103.758 


-.1 


25. 5 


107.872 


107.631 


.2 


26.0 


112.144 


111.523 


.6 


26.5 


1 16.405 


1 15.432 


.8 


27.0 


120.633 


119.355 


1.1 


50.0 


332.800 


299.732 


9.9 


100.0 


682.400 


648.275 


5.0 


150.0 


942.900 


922.095 


2.2 


200.0 


1197.000 


1118.518 


6.6 
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THE LOW DENSITY SAMPLE 
HEAT CAPACITIES IN MJ/OK/GM 



TEMPERATURE 


C ACTUAL 


TARASOV C 


PER CENT 
DIFFERENCE 


2.5 


.288 


.288 


.0 


3.0 


.523 


.520 


.6 


3. 5 


.855 


.855 


-. 1 


4.0 


1.298 


1.304 


-.5 


4.5 


1 .865 


1 .871 


-.3 


5.0 


2.566 


2.564 


• 1 


5. 5 


3.405 


3.386 


.6 


6.0 


4.386 


4.345 


.9 


6.5 


5.510 


5.448 


1 . 1 


7.0 


6.776 


6.697 


1.2 


7. 5 


8. 182 


8.098 


1 .0 


8.0 


9.72 4' 


9.651 


.7 


8.5 


1 1 .397 


1 1 .356 


.4 


9.0 


13.197 


13.210 


-. 1 


9. 5 


13.118 


15.209 


-.6 


10.0 


17.153 


17.549 


-1.1 


10. 5 


19.293 


19.624 


-1.7 


11.0 


21.546 


22.023 


-2.2 


11.5 


23.892 


24.554 


-2.8 


12.0 


26.331 


27.195 


-3.3 


12. 5 


28.869 


29.945 


-3.8 


13.0 


31.473 


32.798 


-4.2 


13. 5 


34.1 69 


35.746 


-4.6 


14.0 


36.945 


38.784 


-5.0 


14. 5 


39.799 


41 .906 


-5.3 


15.0 


42.732 


45. 106 


-5.6 


15. 5 


45.742 


48.578 


-5.8 


16.0 


48.832 


51.716 


-5.9 


16. 5 


52.002 


55.117 


-6.0 


17.0 


55.255 


50.575 


-6.0 


17.5 


58. 592 


62.085 


-6.0 


18.0 


62.017 


65.643 


-5.8 


18. 5 


65.532 


69.246 


-5.7 


19.0 


69. 141 


72.888 


-5.4 


19. 5 


72.845 


76.568 


-5.1 


20.0 


76.647 


80.280 


-4.7 


20.5 


80.549 


84.022 


-4.3 


21.0 


84.551 


87.791 


-3.8 


21.5 


88.651 


91.584 


-3.3 


22.0 


92.348 


95.399 


-2.7 


22. 5 


97.137 


99.233 


-2.2 


23.0 


101.51 1 


105.083 


-1 .S 


23.5 


105.961 


106.948 


-.9 


24.0 


110.476 


110.826 


-.3 


24.5 


115.038 


114.714 


.3 


25.0 


119.628 


118.612 


.8 


25.5 


124.222 


122.517 


1.4 


26.0 


128.791 


126.429 


1.0 


26.5 


133.300 


130.345 


2.2 


27.0 


137.710 


134.265 


2.5 


50.0 


332.800 


310.235 


6.8 


100.0 


682.400 


652.453 


4.4 


150.0 


942.900 


924.277 


2.0 


200.0 


1 197.000 


1 119.869 


6.4 



106 



THE 0 0/0 CRYSTALINE (AMORPHOUS) SAMPLE 



HEAT 


CAPACITIES 


IN MJ/0K/5M 




TEMPERATURE 


C ACTUAL 


TARASOV C 


PER CENT 








DIFFERENCE 


2.5 


.427 


.428 


_ o 
• 


3.0 


.788 


.781 


1 .0 


3. 5 


1.296 


1.295 


. 1 


4.0 


1 .974 


1.983 


-.4 


4.5 


2.84 1 


2.851 


-.3 


5 . 0 


3.908 


3.904 


. 1 


5. 5 


5.179 


5.148 


.6 


6.0 


6.650 


6.590 


1 .0 


6.5 


0.341 


0.235 


1.3 


7.0 


10.222 


1 0.039 


1.3 


7.5 


12.292 


12.152 


1.1 


6.0 


1 4 . 542 


14.421 


.8 


8.5 


16.959 


16.891 


.4 


9.0 


19.531 


1 9.552 


-.1 


9. 5 


22.244 


22.394 


-.7 


10.0 


25.084 


25.403 


-1.3 


10. 5 


20.040 


28.565 


-1 .9 


11.0 


31.097 


31.866 


-2.5 


11.5 


34.246 


35.291 


-3.1 


12.0 


37.475 


50.827 


— 3.6 


12. 5 


40.776 


42.460 


-4.1 


13.0 


44.140 


46.177 


-4.6 


13.5 


47.563 


49.967 


-5.1 


14.0 


■ 51 .041 


53.819 


-5.4 


14. 5 


54.570 


57.723 


-5.8 


15.0 


58. 151 


61.672 


-6.1 


15.5 


61.783 


65.655 


-6.3 


16.0 


65.471 


69.668 


-6.4 


16.5 


69.216 


73.704 


-6.5 


17.0 


73.025 


77.756 


-6.3 


17.5 


76.903 


81.822 


-6.4 


18.0 


80.856 


85.895 


-6.2 


18.5 


84.891 


89.974 


-6.0 


19.0 


89.01 5 


94.053 


- 5.7 


19.5 


95.234 


98.132 


-5.3 


20.0 


97.553 


102.207 


-4.8 


20.5 


101.975 


106.277 


-4.2 


21.0 


106. 502 


110.339 


-3.6 


21.5 


111.153 


1 14.393 


-2.9 


22.0 


115.064 


1 1 8.437 


-2.2 


22.5 


120.685 


122.470 


-1.5 


23.0 


125.584 


126.492 


-.7 


23. 5 


130.541 


130.501 


.0 


24.0 


135.533 


134.498 


.8 


24. 5 


140. 527 


138.481 


1 .5 


25. 0 


145. 481 


142.451 


2.1 


25.5 


150.348 


146.408 


2.6 


26.0 


155.068 


150.351 


3.u 


26. 5 


159.569 


154.280 


3.3 


27.0 


163.770 


158.195 


3.4 


50.0 


332.800 


327.092 


1.7 


100.0 


682.400 


659.158 


3.4 


150.0 


942.900 


927.779 


1 .6 


200.0 


1 197.000 


1 122.037 


6.3 
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